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We study the processes e + e~ — > K + K~ tt + -k~ j, K + K~ 7r°7r°7, and K + K~ K + K~ 7, where the 
photon is radiated from the initial state. About 84000, 8000, and 4200 fully-reconstructed events, 
respectively, are selected from 454 fb _1 of BABAR data. The invariant mass of the hadronic final 
state defines the e + e~ center-of-mass energy, so that the K + K ~7r + 7r~7 data can be compared with 
direct measurements of the e + e~ — > K + K~-k + -k~ reaction. No direct measurements exist for the 
e + e~ — > K + K~ ir°-K° or e + e~ — > K + K~ K + K~ reactions, and we present an update of our previous 
result based on a data sample that is twice as large. Studying the structure of these events, we find 
contributions from a number of intermediate states, and extract their cross sections. In particular, 
we perform a more detailed study of the e + e~ — > 0(lO2O)7T7r7 reaction, and confirm the presence of 
the y(2175) resonance in the 0(1O2O)/ O (98O) and A" + iC"/o(980) modes. In the charmonium region, 
we observe the J/ip in all three final states and in several intermediate states, as well as the tp(2S) 
in some modes, and measure the corresponding products of branching fraction and electron width. 

PACS numbers: 13.66.Bc, 14.40.-n, 13.25.Jx 



I. INTRODUCTION 

Electron-positron annihilation at fixed center-of-mass 
(cm.) energies has long been a mainstay of research in 
elementary particle physics. The idea of utilizing initial- 
state radiation (ISR) to explore e + e~ reactions below 
the nominal cm. energies was outlined in Ref. [l|, and 
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discussed in the context of high-luminosity 4> and B fac- 
tories in Refs. [2H4|. At high cm. energies, e + e _ anni- 
hilation is dominated by quark-level processes producing 
two or more hadronic jets. Low-multiplicity processes 
dominate below or around 2 GeV, and the region near 
the charm threshold, 3.0-4.5 GeV, features a number of 
resonances Thus, studies with ISR events allow us to 
probe a wealth of physics topics, including cross sections, 
spectroscopy and form factors. Charmonium and other 
states with J PC = 1 can be observed, and interme- 
diate states may contribute to the final state hadronic 
system. Measurements of their decay modes and branch- 
ing fractions are important to an understanding of the 
nature of such states. 

Of particular current interest (see Ref. @) is the 
y(2175) state observed to decay to 0(1O2O)/ O (98O) in 
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our previous study Q and confirmed by the BES [|| and 
Belle [9| Collaborations. With twice the integrated lumi- 
nosity (compared to Ref. Q) in the present analysis, we 
perform a more detailed study of this structure. 

The study of e + e~ — > hadrons reactions in data is also 
critical to hadronic-loop corrections to the muon mag- 
netic anomaly, = (g^ — 2)/2. The theoretical predic- 
tions of this anomaly rely on these measurements [Io| . 
Improving this prediction requires not only more precise 
measurements, but also measurements from threshold to 
the highest cm. energy possible. In addition, all the 
important sub-processes should be studied in order to 
properly incorporate possible acceptance effects. Events 
produced via ISR at B factories provide independent and 
contiguous measurements of hadronic cross sections from 
the production threshold to a cm. energy of ~5 GeV. 
With more data we also are able to reduce systematic 
uncertainties in the cross section measurements. 

The cross section for the radiation of a photon of en- 
ergy E 7 in the cm. frame, followed by the production 
of a particular hadronic final state /, is related to the 
corresponding direct e + e~ — > / cross section <x/(s) by 



da lf (so,x) 
dx 



W(s ,x) ■ a f (s (l - x)) 



(1) 



where ^/~sq is the nominal e + e~ cm. energy, x = 2E~ t / \J~sq 
is the fraction of the beam energy carried by the ISR 
photon, and E c m . = \J So(l — x) = yfs is the effective cm. 
energy at which the final state / is produced. The proba- 
bility density function W{s$,x) for ISR photon emission 
has been calculated with better than 1% precision (see, 
e.g. Ref. Q). It falls rapidly as E 1 increases from zero, 
but has a long tail, which in combination with the in- 
creasing oy(so(l — x)) produces a sizable event rate at 
very low E c m .. The angular distribution of the ISR pho- 
ton peaks along the beam directions. For a typical e + e~ 
detector, around 10-15% of the ISR photons fall within 
the experimental acceptance Q • 

Experimentally, the measured invariant mass of the 
hadronic final state defines E c m .. An important feature 
of ISR data is that a wide range of energies is scanned 
continuously in a single experiment, so that no structure 
is missed, and the relative normalization uncertainties in 
data from different experiments are avoided. Further- 
more, for large values of x the hadronic system is col- 
limated, reducing acceptance issues and allowing mea- 
surements down to production threshold. The mass res- 
olution is not as good as the typical beam energy spread 
used in direct measurements, but resolution and absolute 
energy scale can be monitored by means of the measured 
values of the width and mass of well-known resonances, 
such as the J/ip produced in the reaction e + e~ — > J/ipJ- 
Backgrounds from e + e~ — > hadrons events at the nomi- 
nal y^so and from other ISR processes can be suppressed 
by a combination of particle identification and kinematic 
fitting techniques. Studies of e + e~ — > /i + /i~7 and sev- 
eral multi-hadron ISR processes using BABAR data have 
been performed [3, [LlT - [l7j , demonstrating the viability of 



such measurements. These analyses have led to improve- 
ments in background reduction procedures for more rare 
ISR processes. 

The K + K~ir + Tr~ final state has been measured di- 
rectly by the DM1 Collaboration [H| for ^/s <2.2 GeV, 
and we have previously published ISR measurements of 
the K+K-tt+tt- and R + R~ R + R~ final states [lj] for 
E c .m. < 4.5 GeV. Later we reported an updated mea- 
surement of the R + R ~7r + 7r~ final state with a larger 
data sample, together with the first measurement of the 
K + K~Tr°Tr° final state, in which we observed a structure 
near threshold in the </>/o intermediate state 0j- 

In this paper we present a more detailed study of these 
two final states along with an updated measurement of 
the K + K~ K + K~ final state. In all cases we require 
the detection of the ISR photon and perform a set of 
kinematic fits. We are able to suppress backgrounds suf- 
ficiently to study these final states from their respective 
production thresholds up to E c . m . =5 GeV. In addition 
to measuring the overall cross sections, we study the in- 
ternal structure of the final states and measure cross sec- 
tions for a number of intermediate states that contribute 
to them. We also study the charmonium region, measure 
several J/tp and ^{2S) products of branching fraction 
and electron width, and set limits on other states. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage rings at the SLAC National Accelerator Labora- 
tory. The total integrated luminosity used is 454.2 fb~ , 
which includes 413.1 fb~ collected at the T(4S) peak. 
v /s^ = 10.58 GeV, and 41.1 fb _1 collected at about 
v /i^= 10.54 GeV. 

The BABAR detector is described elsewhere [l^. In 
the present work, we use charged-particle tracks recon- 
structed in the tracking system, which is comprised of a 
five double-sided-layer silicon vertex tracker (SVT) and 
a 40- layer drift chamber (DCH) in a 1.5 T axial magnetic 
field. Separation of charged pions, kaons, and protons is 
achieved using a combination of Cherenkov angles mea- 
sured in the detector of internally-reflected Cherenkov 
light (DIRC) and specific-ionization measurements in the 
SVT and DCH. For the present study we use a kaon iden- 
tification algorithm that provides 90-95% efficiency, de- 
pending on momentum, and pion and proton rejection 
factors in the 20-100 range. Photon and electron ener- 
gies are measured in a CsI(Tl) electromagnetic calorime- 
ter (EMC). We use muon identification provided by an 
instrumented flux return (IFR) to select the ^( + /j, _ 7 final 
state used for photon efficiency studies. 

To study the detector acceptance and efficiency, we 
use a simulation package developed for radiative pro- 
cesses. The simulation of hadronic final states, includ- 
ing K+K-tt+tt-j, K+K-n n j and R+R-R+R-j, 
is based on the approach suggested by Czyz and Kiihn 
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[20| . Multiple soft-photon emission from the initial- 
state charged particles is implemented with a structure- 
function technique [U [22[, and photon radiation from 
the final-state particles is simulated by the PHOTOS 
package |23[. T he p recision of the radiative corrections 
is about 1% [H [If. 

We simulate the two K + K~irir (7r + 7r~, 7r°7r n ) final 
states uniformly in phase space, and also according to 
models that include the 0(1020) -t K + K~ and/or 
/o(980) -> 7T7T channels. The K + K~K + K~ final state 
is simulated according to phase space, and also includ- 
ing the 4> — > K + K~ channel. The generated events are 
subjected to a detailed detector simulation [24]], and we 
reconstruct them with the same software chain used for 
the experimental data. Variations in detector and back- 
ground conditions over the course of the experiment are 
taken into account. 

We also generate a large number of potential back- 
ground processes, including the ISR reactions e + e~ — > 
7r + 7r~7r + 7r~7, e + e~ — > 7r + 7r~7r°7r°7, and e + e~ — > 
KsKttj, which can contribute due to particle misidenti- 
fication. We also simulate e + e~ — > ^777, e + e~ — » 07r°7, 
and e + e~ — > 7r + 7r~7r°7, which have larger cross sections 
and can contribute background via missing or spurious 
tracks or photons. In addition, we study non-ISR back- 
grounds resulting from e + e~— > qq (q — u,d,s,c) gener- 
ated using JETSET 12511 and from e + e~^T + r~ gener- 
ated using KORALB [2q. The cross sections for these 
processes are known to about 10% accuracy or better, 
which is sufficiently precise for the purposes of the mea- 
surements in this paper. The contribution from T(AS) 
decays is found to be negligible. 



III. EVENT SELECTION AND KINEMATIC FIT 

In the selection of candidate events, we consider photon 
candidates in the EMC with energy above 0.03 GeV, and 
charged-particle tracks reconstructed in either or both of 
the DCH and SVT, that extrapolate within 0.25 cm of 
the collision axis in the transverse plane and within 3 cm 
of the nominal collision point along this axis. We require 
a photon with cm. energy Ej > 3 GeV in each event, and 
either four charged-particle tracks with zero net charge 
and total momentum roughly (within 0.3 radians) oppo- 
site to the photon direction, or two oppositely-charged 
tracks that combine with other photons to roughly bal- 
ance the high-energy photon momentum. We assume 
that the photon with the largest value of E 1 is the ISR 
photon. We fit the set of charged-particle tracks to a 
common vertex and use this as the point of origin in 
calculating the photon direction(s). If additional well- 
reconstructed tracks exist, the nearest four (two) to the 
interaction region are chosen for the four-track (two- 
track) analysis. Most events contain additional soft pho- 
tons due to machine background or interactions in the 
detector material. 

We subject each candidate event to a set of constrained 



kinematic fits and use the fit results, along with charged- 
particle identification, both to select the final states of in- 
terest and to measure backgrounds from other processes. 
The kinematic fits use the ISR photon direction and en- 
ergy along with the four-momenta and covariance ma- 
trices of the initial e + e~ and the set of selected tracks 
and photons. The ISR photon energy and position are 
additionally aligned and calibrated using the /i + /z~7 ISR 
process, since the two well-identified muons predict pre- 
cisely the position and energy of the photon. This process 
is also used to identify and measure data - Monte Carlo 
(MC) simulation differences in the photon detection ef- 
ficiency and resolution. The fitted three-momentum for 
each charged-particle track and the photon are used in 
further kinematical calculations. 

For the four-track event candidates the fits have four 
constraints (4C). We first fit to the 7r + 7r~7r + 7r~ hypothe- 
sis, obtaining the chi-squared value xln- H the four tracks 
include one identified K + and one identified K~ , we fit 
to the K + K ~tt + it~ hypothesis and retain the event as 
a K + K ~ir + 7r~ candidate. For events with one identi- 
fied kaon, we perform fits with each of the two oppositely 
charged tracks given the kaon hypothesis, and the combi- 
nation with the lower x\k2k ^ s retained if its value is less 
than xI-k- ^ the event contains three or four identified 
K , we fit to the K + K~K + K~ hypothesis and retain 
the event as a K + K~ K + K~ candidate with chi-squared 
value \\ K . 

For the events with two charged-particle tracks and five 
or more photon candidates, we require that both tracks 
be identified as kaons to suppress background from ISR 
7T + 7r~7r 7r° and K^K^if* events. We then pair all non- 
ISR photon candidates and consider combinations with 
invariant mass within ±30 MeV/c 2 of the ir° mass [5j as 
7T° candidates. We perform a six-constraint (6C) fit to 
each set of two non-overlapping ir° candidates, the ISR 
photon, the two charged-particle tracks, and the beam 
particles. Both 7r° candidates are constrained to the 7r° 
mass, and we retain the combination with the lowest chi- 
squared value, X2K2tt°- 

IV. THE K+R-tt+tt- FINAL STATE 
A. Final Selection and Backgrounds 

The X2K2-W distribution in data for the K + K ~Tr + n~ 
candidates is shown in Fig. [Tj (points); the open his- 
togram is the distribution for the simulated K + K ~ir + ir~ 
events. The distributions are broader than those for 
a typical 4C % 2 distribution due to higher order ISR, 
and the experimental distribution has contributions from 
background processes. The simulated distribution is nor- 
malized to the data in the region x\k2-k < - 10 where the 
contributions of the backgrounds and radiative correc- 
tions do not exceed 10%. 

The shaded histogram in Fig. [1] represents the back- 
ground from non-ISR e + e~ — > qq events obtained from 
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FIG. 1: Distribution of x 2 from the four-constraint fit for 
K + K~ 7r + 7r~ candidates in the data (points). The open his- 
togram is the distribution for simulated signal events, nor- 
malized as described in the text. The shaded, cross-hatched, 
and hatched regions represent, respectively, the background 
from non-ISR events, from the ISR KsKn process, and back- 
grounds with dominant contribution from mis-identified ISR 
47r events. Signal and control regions are indicated. 

the JETSET simulation. It is dominated by events with 
a hard tt° that results in a fake ISR photon. These events 
otherwise have kinematics similar to the signal, resulting 
in the peaking structure at low values of x\k2-k ■ We eval- 
uate this background in a number of E c m . ranges by com- 
bining the ISR photon candidate with another photon 
candidate in both data and simulated events, and com- 
paring the 7T° signals in the resulting 77 invariant mass 
distributions. The simulation gives an E c m . -dependence 
consistent with the data, so we normalize it using an over- 
all factor. The cross-hatched region in Fig. Q] represents 
e + e~— > KsKitj events with K$ — > tt + tt~ decays close 
to the interaction region, and one pion mis-identified as 
a kaon. The process has similar kinematics to the signal 
process, and a contribution of about 1% is estimated us- 
ing the cross section measured in our previous study [l6[ ■ 
The hatched region represents the contribution from ISR 
e + e~— > 7r + 7r~7T + 7r~ events with one or two misidentified 
pions; this process contributes mainly at low \ 2 values. 
We estimate the contribution as a function of E c m . from 
a simulation using the cross section value and shape from 
our previous study |13j |. 

All remaining background sources are either negligi- 
ble or give a X 2 k2tt distribution that is nearly uniform 
over the range shown in Fig. [T] We define the signal re- 
gion by requiring x\k2-h < - 30 j an d estimate the sum of 
the remaining backgrounds from the difference between 
the number of data and simulated entries in the control 
region, 30 < x\k2-k < 60, as shown in Fig. Q] The back- 
ground contribution to any distribution other than x 2 is 



estimated as the difference between the distributions in 
the relevant quantity for data and MC events from the 
control region of Fig. [1] normalized to the difference be- 
tween the number of data and MC events in the signal re- 
gion. The non-ISR background is subtracted separately. 
The signal region contains 85598 data and 63784 simu- 
lated events; the control region contains 9684 data and 
4315 simulated events. 
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m(K + K~7tV) (GeV/c 2 ) 

FIG. 2: The invariant mass distribution for K + K~ tt + tv~ 
candidates in the data (points): the shaded, cross-hatched 
and hatched regions show, respectively, the non-ISR back- 
ground from JETSET simulation, the KsKn background 
with a small contribution from the control region of Fig.[TJ and 
the dominant contribution resulting from ISR mis-identified 
7r + 7r _ 7r + 7r _ events. 

Figure [2] shows the K + K~ ~ir + ir invariant mass dis- 
tribution from threshold up to 5.0 GeV/c 2 for events in 
the signal region. Narrow peaks are apparent at the J ftp 
and ip(2S) masses. The shaded histogram represents the 
qq background, which is negligible at low mass but dom- 
inates at higher masses. The cross-hatched region rep- 
resents the background from the KgKn channel (which 
exhibits a 0(1680) peak [lj|) and from the x 2 control 
region. The hatched region represents the contribution 
from mis-identified ISR 7r + 7r~7r + 7r~, and is dominant for 
masses below 3.0 GeV/c 2 . The total background is 6-8% 
at low mass, but accounts for 20-25% of the observed dis- 
tribution near 4 GeV/ c 2 , and increases further for higher 
masses. 

We subtract the sum of backgrounds in each mass in- 
terval to obtain the number of signal events. Considering 
uncertainties in the cross sections for the background pro- 
cesses, the normalization of events in the control region, 
and the simulation statistics, we estimate a systematic 
uncertainty on the signal yield that is 2% or less in the 
1.6-3.3 GeV/c 2 mass region, but increases linearly to 10% 
in the 3.3-5.0 GeV/c 2 region, and is about 20% for the 
masses below 1.6 GeV/c 2 . 
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FIG. 3: (a) The invariant mass distributions for 

K + K~ ir + iv~ MC events that are simulated uniformly in 
phase space, reconstructed in the signal (open) and control 
(hatched) regions of Fig. [1] (b) net reconstruction and selec- 
tion efficiency as a function of mass obtained from this simu- 
lation (the curve represents a third-order polynomial fit). The 
dashed curve is obtained for the (/>(1020)7r + 7r~ final state. 



B. Selection Efficiency 

The selection procedure applied to the data is also 
applied to the simulated signal samples. The resulting 
K + K~~ir + ir~ invariant-mass distributions in the signal 
and control regions are shown in Fig.[3][a) for the uniform 
phase space simulation. This model reproduces the ob- 
served distributions of kaon and pion momenta and polar 
angles. A broad, smooth mass distribution is chosen to 
facilitate the estimation of the efficiency as a function of 
mass. We divide the number of reconstructed simulated 
events in each mass interval by the number generated in 
that interval to obtain the efficiency shown by the points 
in Fig. [3jb). The result of fitting a third-order polyno- 
mial to the points is used for further calculations. We 
simulate events with the ISR photon confined to the an- 
gular range 20-160° with respect to the electron beam 
in the e + e~ cm. frame; this angular range is wider than 
the actual EMC acceptance. The calculated efficiency is 
for this fiducial region, and includes the acceptance for 
the final-state hadrons, the inefficiencies of the detector 
subsystems, and event loss due to additional soft-photon 
emission. 

The simulations including the 0(lO2O)7r + 7r~ and/or 
K + K~ fo(980) channels give very different mass and an- 
gular distributions in the K + K~ir + Tr~ rest frame. How- 
ever, the angular acceptance is quite uniform for ISR 
events (see Ref. [13]), and the efficiencies are within 1% of 
those from the uniform phase space simulation, as shown 
by the dashed curve in Fig. 0HJb) for the 0(lO2O)7r + 7r~ 



final state. 

To study possible mis-modeling of the acceptance, 
we repeat the analysis with tighter requirements. All 
charged tracks are required to lie within the DIRC accep- 
tance, 0.45 < 6* c h < 2.4 radians, and the ISR photon must 
not appear near the edges of the EMC, 0.35 <#isr< 2.4 
radians. The fraction of selected data events satisfying 
the tighter requirements differs from the simulated ratio 
by 1.5%. We take the sum in quadrature of this variation 
and the 1% model variation (2% total) as the systematic 
uncertainty due to acceptance and model dependence. 
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The e + e -^K + K tx + tx cross section as a function 
cm. energy measured with ISR data at BABAR (dots). 



FIG. 4 

of e + e 

The direct measurements from DM1 [18| are shown as the 
open circles. Only statistical errors are shown. 



Our data sample contains about 3000 events in the J/tp 
peak. Comparing this number with and without selection 
on X2K2tr we nn d l ess than a 1% difference between data 
and MC simulation due to mis-modeling of the shape of 
the X2K 2it distribution. This value is taken as an estimate 
of the systematic uncertainty associated with the x\k2tx 
selection criterion. To measure tracking efficiency, we 
consider data and simulated events that contain a high- 
energy photon and exactly three charged-particle tracks, 
which satisfy a set of kinematical criteria, including a 
good x 2 from a kinematic fit to the 7r + 7r~7r + 7r~ hypoth- 
esis, assuming one missing pion track in the event. We 
find that the simulated track-finding efficiency is over- 
estimated by (0.75 ± 0.25)% per track, so we apply a 
correction of +(3 ± 1)% to the signal yield. 

The kaon identification efficiency is studied in 
BABAR using many different test processes (e.g. e + e~— > 
0(1020)7 —> K + K~j), and we conservatively estimate a 
systematic uncertainty of ±1.0% per kaon due to data- 
MC differences in our kaon momentum range. 

The data-MC simulation correction due to ISR photon 
detection efficiency was studied with a sample of e + e~ — > 
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TABLE I: Summary of the cross section measurements for e + e — > K + K -k + tt . Errors are statistical only. 
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the mass spectra, obtained from the MC simulation with 
and without extra-soft-photon (ISR and FSR) radiation, 
we extract i?(E c m ), which gives a correction less than 
1%. Our data, calculated according to Eq.[21 include vac- 
uum polarization (VP) and exclude any radiative effects, 
as is conventional for the reporting of e + e~ cross sections. 
Note that VP should be excluded and FSR included for 
calculations of a^. From data-simulation comparisons 
for the e + e _ — > /i + /.t~7 events we estimate a systematic 
uncertainty on dC of 1% [l7| . 

We show the cross section as a function of E c m . in 
Fig. [4] with statistical errors only in comparison with the 
direct measurements from DM1 [l8j], and list our results 
in Table |U The results are consistent with our previous 
measurements for this reaction 0, [H[ , but have increased 
statistical precision. Our data lie systematically below 
the DM1 data for E c m . above 1.9 GeV. The systematic 
uncertainties, summarized in Table UH affect the normal- 
ization, but have little effect on the energy dependence. 

The cross section rises from threshold to a peak value 



7 events and was found to be +(1.0 ± 0.5)%. 



C. Cross Section for e + e — y K + K 7r + 7r 

We calculate the e + e~^K + K~ir + ir~ cross section as 
a function of the effective cm. energy from 

x _ ^2K27r7(E c .m.) 



rf£(E c . m .) • e2K27r(E c .m.) ' ^(E c .m.) 

(2) 

where E c m . = m,2K27rC 2 with m2K2-n the measured in- 
variant mass of the K + K~ n + ti~ system, dN2K2-n^ the 
number of selected events after background subtraction 
in the interval cffi c . m ., £2x271- (E c . m .) the corrected detec- 
tion efficiency, and R a radiative correction. 

We calculate the differential luminosity cLC(E c . m .) in 
each interval e?E c m , with the photon in the same fiducial 
range as that used for the simulation, using the simple 
leading order (LO) formula described in Ref. 121. From 
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of about 4.6 nb near 1.86 GeV, then generally decreases 
with increasing energy. In addition to narrow peaks at 
the J/tp and ip(2S) mass values, there are several possible 
wider structures in the 1.8-2.8 GeV region. Such struc- 
tures might be due to thresholds for intermediate reso- 
nant states, such as </>/o(980) near 2 GeV. Gaussian fits 
to the distributions of the mass difference between gener- 
ated and reconstructed MC data yield K + K ~n + ir~ mass 
resolution values that vary from 4.2 MeV/c 2 in the 1.5- 
2.5 GeV/c 2 region to 5.5 MeV/c 2 in the 2.5-3.5 GeV/c 2 
region. The resolution functions are not purely Gaus- 
sian due to soft-photon radiation, but less than 10% of 
the signal is outside the 0.025 GeV/c 2 mass interval used 
in Fig. 01 Since the cross section has no sharp structure 
other than the J/ip and ip(2S) peaks discussed in Sec. lIXI 
below, we apply no correction for mass resolution. 



D. Substructures in the K~*~ K tt + -k Final State 

Our previous study 0, E3] showed evidence for many 
intermediate resonances in the K + K~ir + ir~ final state. 
With the larger data sample used here, these can be 
seen more clearly and, in some cases, studied in detail. 
Figure EJa) shows a plot of the invariant mass of the 
K~ir + pair versus that of the K + ir~ pair. Signal for 
the i*sT*(892)° is clearly visible. Figure Mb) shows the 
K^w^ mass distribution (two entries per event) for all 
selected K + K~ir + ir~ events. As we show in our pre- 
vious study JtJ, the signal at about 1400 GeV/c 2 has 
parameters consistent with K^il^SO) ■ Therefore, we 
perform a fit to this distribution using P- and D-wave 
Breit-Wigner (BW) functions for the K*° and K^ sig- 
nals, respectively, and a third-order polynomial function 
for the remainder of the distribution, taking into account 
the Ktt threshold. The fit result is shown by the curves 
in Fig. Mb). The fit yields a K*° signal of 53997 ± 526 
events with m(K*°) = 0.8932 ± 0.0002 GeV/c 2 and 



TABLE II: Summary of corrections and systematic uncer- 
tainties for the e + e~ — ¥ K + K~-k + -k~ cross section measure- 
ments. The total correction is the linear sum of the contribu- 
tions, and the total uncertainty is obtained by summing the 
individual uncertainties in quadrature. 



Source 



Correction 



Uncertainty 



Rad. Corrections 




1% 








Backgrounds 




2% 


5 Ec.m. 


< 3.3 


GeV 






2-10 % 


5 E c . m . 


> 3.3 


GeV 


Model Acceptance 




2% 








X2K2-W Distribution 




1% 








Tracking Efficiency 


+3% 


1% 








Kaon ID Efficiency 




2% 








Photon Efficiency 


+1.0% 


0.5% 








ISR Luminosity 




1% 








Total 


+4.0% 


4% 
4-11% 


j E c . m . 
3 E c . m . 


< 3.3 
> 3.3 


GeV 
GeV 



T(K*°) = 0.0521 ± 0.0007 GeV, and a K*° signal of 
4361 + 235 events with m{K^°) = 1.4274+0.0019 GeV/c 2 
and T(Kf) = 0.0902 ± 0.0056 GeV. These values are 
consistent with current world averages for if* (892)° and 
^2(1430)° [5J , and the fit describes the data well, indi- 
cating that contributions from other resonances decaying 
into i^Ti-T, like X* (1410)° and/or if,*(1430) , are small. 
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FIG. 5: (a) Invariant mass of the K~n + pair versus that of 
the K + n~ pair; (b) The K ± ir^ 1 mass distribution (two entries 
per event) for all selected K + K~ -k + -k~ events: the solid line 
represents a fit including two resonances and a polynomial 
background function, which is shown as the hatched region. 

We combine K*°/K*° candidates within the lines in 
Fig. [SJa) with the remaining pion and kaon to ob- 
tain the if*(892)°7r ± invariant mass distribution shown 
in Fig. Mb), and the AT*(892) 7r ± versus K*(892)°K* 
mass plot in Fig. [6ja). The bulk of Fig. [6ta) shows a 
strong positive correlation, characteristic of K*°Kw fi- 
nal states with no higher resonances. The horizontal 
bands in Fig. M a ) correspond to the peak regions of 
the projection plot of Fig. Mb) and are consistent with 
the contribution from the if 1 (1270) and if 1 (1400) reso- 
nances. There is also an indication of a vertical band in 
Fig. M, a )> perhaps corresponding to a K*(892)°K struc- 
ture at ~1.5 GeV/c 2 . The projection plot of Fig. M, c ) 
for events with m(iC*(892) 7r ± ) > 1.5 GeV/c 2 shows the 
enhancement not consistent with phase space behavior. 

We next suppress the K* (892)° Kir contribution by 
considering only events outside the lines in Fig. [5ja). 
In Fig. Efa) the K ± tt + tt^ invariant mass (two entries 
per event) shows evidence of the iVi (1270) and ifi(1400) 
resonances, both of which decay into Kp(770), although 
the latter decay is very weak [5j. In Fig. Mb) we plot 
the 7r + 7r~ invariant mass for events with m(K ± TT + TT^) > 
1.3 GeV/c 2 . There is a strong p(770) -+ tt+tt^ signal, 
and there are indications of additional structures in the 
/ (980) and / 2 (1270) regions. 

The separation of all these, and any other, intermedi- 
ate states involving relatively broad resonances requires 
a partial wave analysis. This is beyond the scope of this 
paper. Instead we present the cross sections for the sum 
of all states that include K* (892)°, K* (1430)° or p(770) 
signals, and study intermediate states that include a nar- 
row <f> or /o resonance. 
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m(KK + ~) (GeV/c 2 ) 

FIG. 6: (a) Invariant mass of the if*(892)°7r ± system versus 
that of the K*(892)°K^ system; (b) the if* (892)%+ projec- 
tion plot of (a); (c) the if *(892)°if =F projection plot of (a) for 
m(K*(892)°Tr ± ) > 1.5 GeV/c 2 . 



sociated with correlated K* (892)° K* (892)° production 
(about 1% of the total number of if * (892)° events), and 
that 1680±343 events correspond to if *(892)°if 2 *(1430)° 
pairs, compared to 4361 ±235, the total number of events 
with a if|(1430)° in the final state. The distribution 
of the events from the if* (1430)° peak shows a strong 
signal at the if* (892)° mass in Fig. [5£b), which con- 
tains 1648 ± 32 events, in agreement with the number 
of if* (892)°if*, (1430)° pairs obtained above. 




0.75 1 1.25 1.5 
m(K"7t + ) (GeV/c 2 ) 



0.75 1 1.25 1.5 
m(K~7l + ) (GeV/c 2 ) 



FIG. 8: The if ~tt + invariant mass distribution correspond- 
ing to the number of if* (892)° (a) and #1(1430)° (b) events 
obtained from the fits to the K + tt~ invariant mass distribu- 
tion for each interval of K~ir + mass. The curves result from 
the fits described in the text. 



N o 800 
> 

O 600 
§ 400 



I 200 













L (°) j 







1 



1.5 



; 2.5 
m(K + rcV) (GeV/c 2 ) 




1.5 

m(7t + 7C ) (GeV/c 2 ) 



FIG. 7: (a) The invariant mass of the K ± tt + tt~ combinations 
with if* (892)° 'Kir events excluded; (b) the tv + tt~ invariant 
mass for events from (a) with the #1(1270) region suppressed 
by requiring m(K ± TT + n~) > 1.3 GeV/c 2 as shown by vertical 
line in (a). 



E. The e+e" K*(892)°Kn, K% (1430)°Ktt and 
K+K-p(770) Cross Sections 

Signals for the if* (892)° and if* (1430)° are clearly vis- 
ible in the if ± 7r T mass distributions in Fig. EJa,b) . To 
extract the number of events with correlated production 
of if*(892)°if*(892)° and K *(892)°if|(1430)° + c.c, we 
perform the same fit as that shown in Fig.JSJb), but to the 
K + ir~ invariant mass distribution in each 0.04 GeV/c 2 
interval of K~ir + invariant mass. From each fit we 
obtain the number of if* (892)° and if* (1430)° events 
and plot these values as a function of K~tt + mass in 
Fig. [HIa) and Fig. [S^b), respectively. The fit to the data 
of Fig. [S^a) indicates that only 548 ± 263 events are as- 




FIG. 9: The e+e" -s- if* (892)°if "tt + cross section, obtained 
from the if* (892)° signal of Fig.^b). 

We perform a fit similar to that shown in Fig. [SJ^b) 
to the data in intervals of K + K ~tt + h~ invariant mass, 
with the resonance masses and widths fixed to the val- 
ues obtained from the overall fit. Since correlated if* 
production is small, we convert the resulting if* yield 
in each interval into a cross section value for e + e~ — ¥ 
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TABLE III: Summary of the cross section measurements for e + e — > K*°(892)K 7r + . Errors are statistical only. 
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FIG. 10: The Kl (1430)°^"7r + cross section, obtained from 
the if| (1430)° signal of Fig. \Mb). 



K*(892)°K-tt+ or J ftT 2 *(1430)° J ft:-7r+, 1 following the pro- 
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FIG. 11: (a) The 7r + 7r~ mass distribution for all selected 
K + K~tt + it~ events with the <j> and K*° regions excluded: the 
solid curve represents a fit as described in the text, and the 
background contribution is shown separately as the hatched 
region; (b) the e + e~ — > K + K~ p(770) cross section obtained 
from the p signal from the fit in each 0.025 GeV cm. energy 
interval. 



cedure described in Sec. IIV CI These cross section values 
take into account only the Kir decay of the if * (892)° and 
the # 2 *(1430)°. 

Note that the e+e" -> K*{m2)°Kn (iC"|(1430)°iC"7r) 
cross section includes a small contribution from the 
K^(U30)°Kn {K* (892)° Kit) channel, because the 
.&:2( 1 430) o i ; s:*(892) final state has not been taken into 
account. These cross sections are shown in Fig. [9] and 



1 The use of charge conjugate reactions is implied throughout the 



paper 
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Fig. EEQl and the e + e~ -> K*(892)°K-ir+ channel is 
listed in Table IIIII for E c m energies from threshold up 
to 4.0 GeV. At higher energies the signals are small 
and contain an unknown, but possibly large, contribu- 
tion from e + e _ — > qq events. There is a rapid rise from 
threshold to a peak value of about 4 nb at 1.84 GeV 
for the e+e~-^ A'*(892)°A _ 7r + cross section, followed 
by a very rapid decrease with increasing energy. There 
are suggestions of narrow structures in the peak region, 
but the only statistically significant structure is the J/ip 
peak, which is discussed below. There are some struc- 
tures in the e + e~ — > K 2 (1430)°A~7r + cross section, but 
the signal size is too small to make any definite state- 
ment. 

The e + e - -» A*(892)°A"7r+ contribution is a large 
fraction of the total A + A~7r + 7r~ cross section at all 
energies above its threshold, and dominates in the 1.8- 
2.0 GeV region. The K + K~p°(770) intermediate state 
makes up the majority of the remainder of the cross 
section. We exclude a small <fi contribution by requir- 
ing \m{K + K~) — m{4>)\ > 0.01 GeV/c 2 , and suppress 
the large K* (892)° contribution by means of the anti- 
selection ImQA^TrF) - 0.892| > 0.035 GeV/c 2 . Fig- 
ure [Ufa) shows the 7r + 7r~ mass distribution for the re- 
maining events. The combinatorial background is rel- 
atively large, and includes a small contribution from 



/o(980) 



7T 7T 



decays. We fit the p(770) signal with a 



single BW (mass and width are fixed to 0.77 GeV/c 2 
and 0.15 GeV, respectively) and a polynomial back- 
ground (contribution shown by the hatched area) in each 
0.025 GeV cm. energy interval. The cross section ob- 
tained is shown in Fig. Illf b). and has no significant 
structures except the J/tjj signal. The uncertainty in the 
/)(770) shape, and also in the background shape, pro- 
vides the largest contribution to the systematic error, 
estimated to be 20-30%. A small contribution to the 
background from /o(980) — > n + ir~ is ignored in the fit, 
which does not result in a significant uncertainty. 



F. The </)(1020)7r+7r~ Intermediate State 

Intermediate states containing narrow resonances can 
be studied more easily. For the E c . m . energy range below 
3.0 GeV, Fig. IT2T a) shows a plot of the invariant mass 
of the 7r + 7r~ pair versus that of the K + K~ pair. Hor- 
izontal and vertical bands corresponding to the p°(770) 
and 4>, respectively, are visible, and there is a concentra- 
tion of entries in the <j) band corresponding to the corre- 
lated production of <fi and /o(980), as demonstrated by 
the open histogram of Fig. [T27b). The (j> signal is clearly 
visible in the K + K~ mass projection of Fig. IT2T c). The 
large contribution from the p(770) is nearly uniform in 
K + K~ mass, and the cross-hatched histogram shows the 
non-A + K~tt + tt~ background estimated from the control 
region in x\k2it- The cross-hatched histogram also shows 
a (j> peak, but this is a small fraction of the events. When 
we subtract this background and fit the remaining data 
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FIG. 12: (a) m(-K + -K~) versus m(K + K~) for all selected 
K + K~ 7r + 7r~ events; (b) the -k + it~ invariant mass projec- 
tions for events in the (j> peak (open histogram), sidebands 
(hatched), and background control region (cross-hatched); (c) 
the K + K~ mass projections for all events (open) and control 
region (cross-hatched); (d) the difference between the open 
histogram and the sum of the other contributions to (b). 



with a double-Gaussian function for the <f> signal, and a 
first-order polynomial function for the non-0 background 
(with a cut-off at the KK threshold), we obtain 3951 ±91 
events corresponding to the 07r + 7r~ intermediate state. 

To study the 07r + 7r~ channel, we select candidate 
events with a K + K~ invariant mass within 10 MeV/c 2 
of the <f> mass, indicated by the inner vertical lines in 
Figs. [T2T a,c). and estimate the non-0 contribution from 
the mass sidebands between the inner and outer vertical 
lines. In Fig. [T2T b') we show the 7r + 7r~ invariant mass dis- 
tributions for <j) candidate events, sideband events, and 
X 2 control region events as the open, hatched and cross- 
hatched histograms, respectively, and in Fig. IT2T d) we 
show the 7r + 7r~ distribution after subtracting the non- 
4> background contributions. We observe a clear, nar- 
row peak in the /o(980) mass region, together with a 
broad enhancement that reaches a maximum at about 
0.6 GeV/c 2 , which could indicate / (600) production. We 
defer a detailed analysis of this distribution to Secs. lIV Gl 
IVTTl and ivnTl 

We obtain the number of e + e~ — > 07t + 7t~ events in 
0.025 GeV/c 2 intervals of the 07r + 7r _ invariant mass by 
fitting the K + K~ invariant mass projection in that in- 
terval after subtracting non-A + A~7r + 7r~ background. 
Each projection is a subset of Fig. [T27 c), where the curve 
represents the fit to the full sample. In each mass inter- 
val, all parameters other than number of events in the 
<f> peak and the normalization of the background distri- 
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TABLE IV: Summary of the cross section measurements for e + e — > <j!>(1020)7r + 7r . Errors are statistical only. 



Earn. (GeV) 


G 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


E c . m . (GeV) 


a 


(nb) 


E c . m . (GeV) 


a (nb) 


1.4875 


0.04 


± 0.01 


1.8375 


0.29 


± 0.07 


2.1875 


0.33 


± 0.04 


2.5375 


0.07 ± 0.02 


1.5125 


0.08 


± 0.03 


1.8625 


0.36 


± 0.07 


2.2125 


0.21 


± 0.04 


2.5625 


0.05 ± 0.02 


1.5375 


0.10 


± 0.03 


1.8875 


0.34 


± 0.06 


2.2375 


0.19 


± 0.03 


2.5875 


0.07 ± 0.02 


1.5625 


0.12 


± 0.03 


1.9125 


0.29 


± 0.06 


2.2625 


0.10 


± 0.04 


2.6125 


0.07 ± 0.02 


1.5875 


0.21 


± 0.04 


1.9375 


0.32 


± 0.06 


2.2875 


0.13 


± 0.03 


2.6375 


0.10 ± 0.02 


1.6125 


0.33 


± 0.05 


1.9625 


0.28 


± 0.05 


2.3125 


0.13 


± 0.03 


2.6625 


0.07 ± 0.02 


1.6375 


0.48 


± 0.06 


1.9875 


0.32 


± 0.05 


2.3375 


0.15 


± 0.03 


2.6875 


0.04 ± 0.01 


1.6625 


0.49 


± 0.06 


2.0125 


0.37 


± 0.05 


2.3625 


0.12 


± 0.03 


2.7125 


0.05 ± 0.01 


1.6875 


0.54 


± 0.07 


2.0375 


0.31 


± 0.05 


2.3875 


0.10 


± 0.03 


2.7375 


0.06 ± 0.01 


1.7125 


0.53 


± 0.07 


2.0625 


0.39 


± 0.05 


2.4125 


0.12 


± 0.02 


2.7625 


0.03 ± 0.01 


1.7375 


0.48 


± 0.07 


2.0875 


0.32 


± 0.05 


2.4375 


0.12 


± 0.02 


2.7875 


0.02 ± 0.01 


1.7625 


0.61 


± 0.08 


2.1125 


0.51 


± 0.05 


2.4625 


0.10 


± 0.02 


2.8125 


0.03 ± 0.01 


1.7875 


0.48 


± 0.07 


2.1375 


0.41 


± 0.05 


2.4875 


0.08 


± 0.02 


2.8375 


0.04 ± 0.01 


1.8125 


0.34 


± 0.07 


2.1625 


0.43 


± 0.05 


2.5125 


0.08 


± 0.02 


2.8625 


0.03 ± 0.01 
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FIG. 13: The e + e~ — > 4>-k + -k~ cross section as a function 
of e + e~ cm. energy obtained by BABAR (dots) and Belle 
(squares) [§]. 



bution are fixed to the values obtained from the overall 
fit. As a check, we also describe the background as a 
linear function, with all parameters free in each mass in- 
terval; the alternative fit yields consistent results with the 
nominal fit to within 5%, which is taken as a systematic 
uncertainty. 

The reconstruction efficiency may depend on the de- 
tails of the production mechanism. Using the two-pion 
mass distribution in Fig. [l"2T d) as input, we simulate the 
7r + 7r~ system as an S-wave composition of two structures 
both described by the BW amplitudes, with parameters 
set to the values obtained in Sec. IVIII The BW ampli- 
tudes represent the / (980) and the bump at 0.(3 GeV/c 2 , 
which we call / (600) (see Sec. IVIII) . We describe the 
(J)tt + it~ mass distribution using a simple model with one 
resonance of mass 1.68 GeV/c 2 and width 0.3 GeV, which 



decays to </>7r + 7r~ or 0/ o (98O) when phase space allows. 
The reconstructed spectrum that results then has a sharp 
increase at about 2 GeV/c 2 due to the c/>/ (980) thresh- 
old. 

We obtain the efficiency as a function of cf>ir + ir~ mass 
by dividing the number of reconstructed events in each 
interval by the number generated; the result is shown in 
Fig. |31 by the dashed curve. Comparison with the solid 
curve in the same figure shows that the model depen- 
dence is weak, giving confidence in the efficiency calcu- 
lation. We calculate the e + e~— ^ <f>ir + ir~ cross section as 
described in Sec. I IV CI and divide by the cf) — > K + K~ 
branching fraction (0.489 [5J). We show our results as 
a function of cm. energy in Fig. 1131 and list them in 
Table IIVI The cross section has a peak value of about 
0.6 nb at about 1.7 GeV, then decreases with increas- 
ing energy until the 0(lO2O)/o(98O) threshold, around 
2.0 GeV. From this point it rises, falls sharply at about 
2.2 GeV, and then decreases slowly. Except in the char- 
monium region, the results at energies above 3 GeV are 
not meaningful due to small signals and potentially large 
backgrounds, and are omitted from Table ITVl Figure [T3l 
displays the cross section up to 4.0 GeV in order to 
show the J/i/j and signals, which are discussed 

in Sec. US 

The cross section obtained is in agreement with our 
previous measurement Q- The cross section measured 
by the Belle Collaboration Q, also shown in Fig. IT31 
presents very similar features, and a general consistency 
with our data, although a small systematic difference at 
higher cm. energies is visible. 

We perform a study of the angular distributions in the 
0(lO2O)7r + 7r _ final state by considering all K + K~n + ir~ 
candidate events with mass below 3 GeV/c 2 in intervals 
of the cosine of each angle defined below, and fitting the 
background-subtracted K + K~ mass projection in each 
interval. The efficiency is nearly uniform in the cosine 
of each angle, and so we study the number of events in 
each interval. We define the (j> production angle, 0^, as 
the angle between the <fi direction and the ISR photon 
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TABLE V: Summary of the e + e~ -> ^(1020)7T7t cross section, dominated by 0(1020) f (980), /o(980) ->■ tvtt, obtained from 
0(lO2O)7r + 7r - events with 0.85 < m(7r + 7r" '-) < 1.1 GeV/c 2 . Err ors are statistical only. 



E c . m . (GeV) 


a (nb) 


E c . m . (GeV) 


a (nb) 


E c . m . (GeV) 


cr (nb) 


E c . m . (GeV) 


a (nb) 


1.8875 


0.00 ± 0.01 


2.1625 


0.54 ± 0.06 


2.4375 


0.11 ± 0.02 


2.7125 


0.04 ± 0.03 


1.9125 


0.01 ± 0.02 


2.1875 


0.38 ± 0.05 


2.4625 


0.11 ± 0.03 


2.7375 


0.04 ± 0.02 


1.9375 


0.16 ± 0.04 


2.2125 


0.19 ± 0.04 


2.4875 


0.08 ± 0.02 


2.7625 


0.03 ± 0.02 


1.9625 


0.15 ± 0.04 


2.2375 


0.19 ± 0.04 


2.5125 


0.07 ± 0.02 


2.7875 


0.03 ± 0.02 


1.9875 


0.19 ± 0.04 


2.2625 


0.10 ± 0.04 


2.5375 


0.06 ± 0.02 


2.8125 


0.02 ± 0.02 


2.0125 


0.32 ± 0.05 


2.2875 


0.15 ± 0.03 


2.5625 


0.05 ± 0.02 


2.8375 


0.05 ± 0.02 


2.0375 


0.28 ± 0.05 


2.3125 


0.14 ± 0.03 


2.5875 


0.07 ± 0.02 


2.8625 


0.03 ± 0.02 


2.0625 


0.38 ± 0.06 


2.3375 


0.16 ± 0.03 


2.6125 


0.07 ± 0.02 


2.8875 


0.02 ± 0.02 


2.0875 


0.35 ± 0.05 


2.3625 


0.14 ± 0.03 


2.6375 


0.07 ± 0.02 


2.9125 


0.04 ± 0.02 


2.1125 


0.60 ± 0.06 


2.3875 


0.07 ± 0.03 


2.6625 


0.07 ± 0.02 


2.9375 


0.01 ± 0.02 


2.1375 


0.50 ± 0.07 


2.4125 


0.11 ± 0.03 


2.6875 


0.03 ± 0.02 


2.9625 


0.01 ± 0.01 





FIG. 14: Distributions of the cosine of (a) the <f) production angle, (b) the pion helicity angle, and (c) the kaon helicity angle 
(see text) for e + e~ — > 07r + 7r~ events: the curves (normalized to the data) represent the distributions expected if the tt + tv~ 
system recoiling against the vector 4> meson is an S-wave system produced in an S-wave orbital angular momentum state. 
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E c . m .(GeV) E cm (GeV) 

FIG. 15: The e + e~— > 4>-k + -k~ cross section derived from the 
K + K~-k + -k~ final state as a function of cm. energy, for (a) 
the 0.85 < m(7r + 7r - ) < 1.1 GeV/c 2 region, dominated by the 
0(1O2O)/ O (98O), and (b) m(ir+-K-) < 0.85 GeV/c 2 . 



direction in the rest frame of the <p-K + iT~ system (i.e., the 
effective e + e~ collision axis). The distribution of cosO^, 
shown in Fig. I14f a). is consistent with the uniform dis- 
tribution expected if the quasi- two-body final state 4>X, 
X — >7r + 7r~, is produced in an S-wave angular-momentum 
state. We define the pion helicity angle, , as that be- 
tween the 7r + and the recoil <f> direction in the tt + tt~ rest 



frame. The kaon helicity angle, Qk+ is defined as that 
between the K + direction and the ISR photon direction 
in the <f> rest frame. The distributions of cos 0^+ and 
cos9 A '+, shown in Figs. IT4Tb) and IT4Tc). respectively, 
are consistent with those expected from scalar (uniform) 
and vector (cos 2 0^+) meson decays, where for the latter 
the <f> retains the helicity of the virtual photon to which 
the 4>X system couples. 



G. The 4>( 1020) / (980) and </>(1020)/ (600) 
Intermediate States 

The narrow /o(980) peak seen in Fig. [T27d) allows the 
selection of a fairly clean sample of </>/o(980) events. We 
repeat the analysis just described with the additional re- 
quirement that the Tr + n~ invariant mass be in the range 
0.85-1.10 GeV/c 2 . A fit to the K+R- mass spectrum 
for this sample, analogous to that shown in Fig. [127c). 
yields about 1350 events; all of these contain a true 
4>, with a small fraction of events with the pion pair 
not produced through the /o(980), but the latter con- 
tribution is relatively small (see discussion in Sec. IVIII) . 
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TABLE VI: Summary of the cross section measurements, dominated by e + e —5- (?i(1020)/o(600), /o(600) — > tt-k process, 



obtained for 


m(7r + 7r") < 0.85 GeV/c 2 . Errors 


are statistical only. 














Ec.m. (GeV) 


a (nb) 


Ec.m. (GeV) 


(7 


(nb) E c 


.m. (GeV) 


a 


(nb) 


E c . m . (GeV) 


a 


(nb) 


1.2875 


0.00 ± 0.01 


1.7125 


0.79 


± 0.11 


2.1375 


0.10 


± 0.04 


2.5625 


0.00 


± 0.01 


1.3125 


0.01 ± 0.01 


1.7375 


0.72 


± 0.10 


2.1625 


0.10 


± 0.04 


2.5875 


0.02 


± 0.01 


1.3375 


0.00 ± 0.01 


1.7625 


0.91 


± 0.12 


2.1875 


0.05 


± 0.03 


2.6125 


0.03 


± 0.01 


1.3625 


0.01 ± 0.01 


1.7875 


0.72 


± 0.11 


2.2125 


0.05 


± 0.03 


2.6375 


0.03 


± 0.02 


1.3875 


0.01 ± 0.01 


1.8125 


0.51 


± 0.10 


2.2375 


0.06 


± 0.03 


2.6625 


0.01 


± 0.01 


1.4125 


0.00 ± 0.01 


1.8375 


0.43 


± 0.10 


2.2625 


0.04 


± 0.02 


2.6875 


0.02 


± 0.02 


1.4375 


0.02 ± 0.01 


1.8625 


0.54 


± 0.11 


2.2875 


0.03 


± 0.02 


2.7125 


0.02 


± 0.02 


1.4625 


0.05 ± 0.02 


1.8875 


0.50 


± 0.09 


2.3125 


0.03 


± 0.02 


2.7375 


0.03 


± 0.03 


1.4875 


0.06 ± 0.02 


1.9125 


0.40 


± 0.09 


2.3375 


0.08 


± 0.02 


2.7625 


0.01 


± 0.02 


1.5125 


0.12 ± 0.04 


1.9375 


0.32 


± 0.08 


2.3625 


0.04 


± 0.02 


2.7875 


0.00 


± 0.01 


1.5375 


0.15 ± 0.04 


1.9625 


0.26 


± 0.07 


2.3875 


0.06 


± 0.02 


2.8125 


0.01 


± 0.02 


1.5625 


0.18 ± 0.04 


1.9875 


0.27 


± 0.07 


2.4125 


0.05 


± 0.02 


2.8375 


0.01 


± 0.02 


1.5875 


0.31 ± 0.06 


2.0125 


0.25 


± 0.06 


2.4375 


0.04 


± 0.02 


2.8625 


0.03 


± 0.02 


1.6125 


0.48 ± 0.08 


2.0375 


0.18 


± 0.05 


2.4625 


0.03 


± 0.01 


2.8875 


0.01 


± 0.02 


1.6375 


0.70 ± 0.09 


2.0625 


0.25 


± 0.05 


2.4875 


0.01 


± 0.01 


2.9125 


0.02 


± 0.02 


1.6625 


0.72 ± 0.09 


2.0875 


0.15 


± 0.05 


2.5125 


0.02 


± 0.01 


2.9375 


0.00 


± 0.01 


1.6875 


0.80 ± 0.10 


2.1125 


0.18 


± 0.05 


2.5375 


0.03 


± 0.01 


2.9625 


0.00 


± 0.01 



By selecting events with the 7r + 7r~ invariant mass be- 
low 0.85 GeV/c 2 , we similarly obtain a sample composed 
mostly of (/)/ (600) events. 

We convert the above two samples of /o(980) and 
/o(600) events in each mass interval into measurements 
of the e+e-^0(lO2O)/ o (98O) and e+e"-)- 0(1020) / (600) 
cross sections as described above, dividing by the /q— > 
ir + ir~ branching fraction of 2/3 to account for /o — > 7r°7r° 
decays. The cross sections are shown in Fig. [ll)]as func- 
tions of cm. energy and are listed in Table |V] and Ta- 
ble EH The 0(1O2O)/ O (98O) cross section behavior near 
threshold does not appear to be smooth, but is more 
consistent with a steep rise to a value of about 0.3 nb 
at 2.0 GeV followed by a slow decrease that is inter- 
rupted by a structure around 2.175 GeV. In contrast, 
the </>(1020)/o(600) cross section has a smooth threshold 
increase to about 0.8 nb, followed by a smooth decrease 
thereafter, and can be interpreted as the (£(1680) reso- 
nance. It is important to note that all structures above 
2.0 GeV seen in Fig. [13] relate only to the / (980) res- 
onance. Possible interpretations of these structures are 
discussed in Sec. lVIIll Again, the cross section values are 
not meaningful for cm. energy above about 3 GeV, ex- 
cept for the J/ip and '0(25') signals, discussed in Sec. IIXI 

V. THE K+K--k°-k° FINAL STATE 
A. Final Selection and Backgrounds 

The K + K~ir°ir° sample contains background from 
the ISR processes e + e~ — > K + K~tt°j and K + K~rjj, 
in which two soft photon candidates from machine- or 
detector-related backgrounds combine with the relatively 
energetic photons from the 7r° or r\ to form two fake it 
candidates. We reduce this background using the angle 



between each reconstructed ir° direction and the direc- 
tion of its higher-energy photon daughter calculated in 
the 7T° rest frame. If the cosines of both angles are larger 
than 0.85, we remove the event. 

Figure [TH] shows the distribution of X2K2tt° f° r 
the remaining candidates together with the simulated 
K + K~ir°Tr° events. Again, the distributions are broader 
than those for a typical 6C % 2 distribution due to higher 
order ISR, and we normalize the histogram to the data 
in the region X 2 K27t<' < 15. The cross-hatched histogram 
in Fig. \W\ represents background from e + e~^rqq events, 
evaluated in the same way as for the K + K~ix + ir~ fi- 
nal state. The hatched region represents the ISR back- 
grounds from final states with similar kinematics. The 
first of these is 7r + 7r~7r°7r°, which yields events with both 
charged pions misidentified as kaons, and the second is 
the KgKir, which yields K$ — > 7r°7r° and a misidentified 
pion. Each contribution is small. 

The dominant background in this case is from resid- 
ual ISR K+R-tt and K+K^-q events, as well as ISR- 
produced K + K~ tt tt°it events. Their net simulated con- 
tribution, indicated by the dashed contour in Fig. I16[ 
is consistent with the data in the high x\k2-k° re gi° n - 
All other backgrounds arc cither negligible or distributed 
uniformly in x\ K27T o- We define the signal region by 
X2K2tt° < 50, which contains 7967 data and 7402 sim- 
ulated events, and a control region by 50 < x\k ! >-k° < 
100, which contains 2007 data and 704 simulated signal 
events. 

Figure [T7] shows the K + K~ir°ir° invariant mass dis- 
tribution from threshold up to 5 GeV/c 2 for events in 
the signal region. The qq background (cross-hatched his- 
togram) is negligible at low masses but yields a signifi- 
cant fraction of the selected events above about 4 GeV/ c 2 . 
The ISR 7T + 7r~7r 7r° contribution (hatched region) is neg- 
ligible except in the 1.5-2.5 GeV/c 2 region. The sum of 
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FIG. 16: Distribution of x 2 from the six-constraint fits 
to K + K~ tt°ty° candidates in the data (points). The open 
histogram is the distribution for simulated signal events, 
normalized as described in the text. The cross-hatched, 
hatched, and dashed regions represent, respectively, the back- 
grounds from non-ISR qq events, ISR-produced 7r + 7r~7r°7r° 
and KsKtv events, and ISR-produced K + K~tv°, K + K~r) and 
K+K -ttW events. 

all other backgrounds, estimated from the control region, 
is the dominant contribution below 2.5 GeV/c 2 and is 
non-negligible everywhere. The total background varies 
from 100% below 1.6 GeV/c 2 to 25% at higher masses. 

We subtract the sum of the estimated background con- 
tributions from the number of selected events in each 
mass interval to obtain the number of signal events. Con- 
sidering uncertainties in the cross sections for the back- 
ground processes, the normalization of events in the con- 
trol region and the simulation statistics, we estimate a 
systematic uncertainty on the signal yield after back- 
ground subtraction of about 5% in the 1.6-3.0 GeV/c 2 
region; this increases linearly from 5% to 15% in the re- 
gion above 3 GeV/c 2 . 

B. Selection Efficiency 

The detection efficiency is determined in the same 
manner as in Sec. IIVBI Figure HUJa) shows the sim- 
ulated K + K~ir°ir° invariant mass distributions in the 
signal and control regions obtained from the phase space 
model. We divide the number of reconstructed events in 
each 0.04 GeV/c 2 mass interval by the number generated 
in that interval to obtain the efficiency estimate shown 
by the points in Fig. [PBTb); a third-order-polynomial fit 
to the efficiency is used in calculating the cross section. 
Again, the simulation of the ISR photon covers a lim- 
ited angular range, which is about 30% wider than the 
EMC acceptance. Simulations assuming dominance of 




m(K + K"7tV) (GeV/c 2 ) 

FIG. 17: Invariant mass distribution for K + K~ tt°tv° can- 
didates in the signal region for data (points). The cross- 
hatched, hatched, and open regions represent, respectively, 
the non-ISR qq background, the contribution from ISR- 
produced 7r + 7T~7r 7r° and KsKir events, and the contribution 
from the other ISR processes described in the text. 

the <j> — » K + K~ and/or the fa— >ir°ir° channels give re- 
sults consistent with those of Fig. H5fb). and we apply 
a 3% systematic uncertainty for possible model depen- 
dence, as in Sec. IIVBI 

We correct for mis-modeling of the track-finding and 
kaon identification efficiencies as in Sec. IIVBI (correc- 
tions of (+1.9 ± 0.6)% and (0 ± 2.0)%, respectively). 
We do not observe any large discrepancy in the shape 
of the X2R"2tt° distribution, and so apply no correction 
for the X2K2tt° < ^0 selection, but introduce 3% as an 
associated systematic uncertainty. We correct the ir°- 
finding efficiency using the procedure described in detail 
inRef. From ISR e + e~->u;7r 7^7r + 7r~7r 7r 7 events 
selected with and without the tt° from the uj decay, we 
find that the simulated efficiency for one ir° is too large by 
(3.0±1.0)%, and we apply a correction of (+6.0 ± 2.0)% 
because of the two 7r°s in each event. 



C. Cross Section for e+e" -»■ K+K--k°-k 

We calculate the cross section for e + e~ K + K~tt°tt° 
in 0.04 GeV E cm . intervals from the analog of Eq.@, us- 
ing the invariant mass of the K + K~ir°ir° system to de- 
termine the cm. energy. We show the results in Fig. Q1J] 
and list the values and statistical errors in Table IVlTI The 
cross section rises to a peak value near 0.8 nb at 2 GeV 
then shows a rapid decrease, which is interrupted by a 
large J/ip signal; the charmonium region is discussed in 
Sec. IIXI below. The drop at 2.2 GeV is similar to that 
seen for the K + K~iy + tt~ final state. Again, the differ- 
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TABLE VII: Summary of the cross section measurements for e + e — > K + K iy°tv°. Errors are statistical only. 
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FIG. 18: (a) Invariant mass distribution for simulated 
K + K~Tv°iT events in the signal (open) and control (hatched) 
regions (see Fig. I16[l ; (b) net reconstruction and selection ef- 
ficiency as a function of mass obtained from this simulation 
(the curve represents the result of a third-order-polynomial 
fit). 



ential luminosity includes corrections for vacuum polar- 
ization that should be omitted for calculations of a^. 

The simulated K + K~ir°ir° invariant mass resolution 
is 8.8 MeV/c 2 in the 1.5-2.5 GeV/c 2 mass range, and 
increases with mass to 11.2 MeV/c 2 in the 2.5-3.5 GeV/c 2 
range. Since less than 20% of the events in a 0.04 GeV 
interval are reconstructed outside that interval, and the 
cross section has no sharp structure other than the J/ip 
peak, we again make no correction for resolution. The 
point-to-point systematic uncertainties are much smaller 
than the statistical uncertainties, and the errors on the 



TABLE VIII: Summary of corrections and systematic un- 
certainties for the e + e~— > K + K^tt°ty° cross section measure- 
ments. The total correction is the linear sum of the contribu- 
tions, and the total uncertainty is obtained by summing the 
individual contributions in quadrature. 



Source 


Correction 


Uncertainty 




Rad. Corrections 




1% 




Backgrounds 




5% , E c . m . < 3 


GeV 






5-15% , Ec.m. > 3 


GeV 


Model Dependence 




3% 




xIk2ttO Distribution 




3% 




Tracking Efficiency 


+1.9% 


0.6% 




Kaon ID Efficiency 




2% 




7T° Efficiency 


+6% 


2% 




ISR-photon Efficiency 


+1.0% 


0.5% 




ISR Luminosity 




1% 





Total 


+8.9% 


7% 


3 E c . m . 


< 3 


GeV 






7-16% 


) E c .m. 


> 3 


GeV 



normalization are summarized in Table IVIIIi along with 
the corrections that were applied to the measurements. 
The total correction is +8.9%, and the total systematic 
uncertainty is 7% at low mass, increasing linearly from 
7% to 16% above 3 GeV/c 2 . 



D. Substructure in the K+ K'tt ^ Final State 

A plot of the invariant mass of the K~tt° pair ver- 
sus that of the K + ir° pair is shown in Fig. l20T a) (two 
entries per event) for the \ 2 signal region after remov- 
ing the ^(1020) contribution by \m{K + K~) - m{4>)\ > 
0.01 GeV/c 2 . Horizontal and vertical bands correspond- 
ing to the if *(892)~ and if *(892)+, respectively, are vis- 
ible. Figure [2UT b) shows as points the sum of the two 
projections of Fig. l2"0Ta); a large if*(892) ± signal is evi- 
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3.5 4 
E c„ (GeV) 

FIG. 19: The e + e~— > K + K~ir°-K° cross section as a function 
of e + e~ cm. energy measured with ISR data at BABAR. The 
errors are statistical only. 




12 12 
m(K + 7t°) (GeV/c 2 ) m(K + 7t°, K"jc°) (GeV/c 2 ) 
FIG. 20: (a) Invariant mass of the K~ir° pair versus that of 
the K + ir° pair in selected if + K~-k°ty° events (two entries per 
event); (b) sum of the projections of (a) (dots, four entries per 
event). The curves represent the result of the fit described in 
the text. 



dent. Fitting this distribution with the function used in 
Sec. IIVE1 we obtain the number of events corresponding 
to K*(892)± (7734 ± 320) and 7^*(1430) ± (793 ± 137) 
production. The K*(U30) ± :K*(892) ± ratio is consis- 
tent with that obtained for neutral K* production in the 
K + K~tt + tt~ channel, but the number of K*(892) ± com- 
binations in the peak is larger than the total number of 
K + K~ir°ir° events (5522). This indicates the presence 
of some number of correlated K* (892) + K* (892)~ pairs. 
Fitting the K~ir° mass distribution in each 0.04 GeV/c 2 
bin of K + ir° invariant mass, we obtain the number of 
2<r*(892) - and K*(U30)~ events shown in Fig. |2"TIa,b), 
The correlated production of *f*(892)+ir*(892)- and 
K*(892)+K*(U30)- is clearly seen, and the fits yield 
1750 ± 60 and 140 ± 49 events, respectively. Note 
that K* (892)+ K* (892)~ accounts for about 30% of all 
K + K~ir°ir° events, in contrast with the K + K~ tt + it~ 
channel, where only 548 ± 263 events (less than 1% of 




0.75 1 1.25 1.5 
m(K + 7t°) (GeV/c 2 ) 

FIG. 21: The number of 7T(892)- (a) and A" 2 *(1430)- (b) 
events obtained from the fits to the K~ir invariant mass 
distributions for each 0.04 GeV/c 2 interval of K + tv° mass. 
The curves result from the fits described in the text. 



0.75 1 1.25 1.5 
m(K + Ji°) (GeV/c 2 ) 
The number of if* (892) ~ 



the total) are found to result from the K*(892)°K*(892) a 
pair production. 

We find no evidence for resonance production in the 
K + K~tt° or K ± 7t tt q subsystems. Since the statistics 
are low in any given mass interval, we do not attempt 
to extract a separate K* (892) + K~ ir° + c.c. cross section. 
The total K + K~ir°ir° cross section is roughly a factor 
of four lower than the K*(892)°K~ir + cross section ob- 
served in the K + K~ir + ir~ final state. This is consistent 
with what is expected from isospin considerations and the 
charged versus neutral K* branching fractions involving 
charged kaons. 



E. The </>(1020)7r°7r° Intermediate State 

The selection of events containing 0(1020) — > K + K~ 
decays follows that in Sec. HVFl Figure l2"27 a) shows the 
plot of the invariant mass of the 7r°7r° pair versus that 
of the K + K~ pair. The 4> resonance is visible as a verti- 
cal band, whose intensity decreases with increasing 7r 7r° 
mass except for an enhancement in the /o(980) region 
(Fig. 122Tb)). The signal is also visible in the K+R- 
invariant mass projection for events in the control re- 
gion, shown in Fig. l22f c). The relative non-0 background 
is smaller than in the K + K~tt + tt~ mode, but there is 
a large background from ISR 07r°, <\>r\ and/or 07r O 7r°7r° 
events, as indicated by the control region histogram 
(hatched) in Fig. I2"27c). The contributions from non- 
ISR and ISR 7r + 7r~7r°7r° events are negligible. Selecting 
4> candidate and side band events as for the K + K~ir + 7r~ 
mode (vertical lines in Figs. I2"27 a.c)). we obtain the 7r°7r° 
mass projections shown as the open and cross-hatched 
histograms, respectively, in Fig. [22T b). Control region 
events (hatched histogram) are concentrated at low mass 
values in Fig. I2"27b). and a peak corresponding to the 
/o(980) is visible over a relatively low background. 

In Fig. I2"27d) we show the 7r°7r° mass distribution as- 
sociated with (j) production after subtraction of all back- 
ground contributions. The distribution is consistent in 
shape with that of Fig. H"2T d). but with a data sample 
which is about six times smaller. 
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TABLE IX: Summary of the e + e~ -> 4>(W2Q)ttit cross section, dominated by 0(1O2O)/ O (98O), /o(980) -> tttt, obtained from 
(?!>(1020)7r 7r events with 0.85 < m(7r°7r°) < 1.1 GeV/c 2 . Errors are statistical only. 



Earn. (GeV) 


a (nb) 


E„n. (GeV) 


a (nb) 


E c . m . (GeV) 


a (nb) 


E c . m . (GeV) 


a (nb) 


1.9000 


0.15 ± 0.07 


2.1000 


0.45 ± 0.11 


2.3200 


0.12 ± 0.06 


2.9200 


0.02 ± 0.01 


1.9400 


0.14 ± 0.06 


2.1400 
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2.4000 
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3.0800 
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2.1800 
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2.4800 
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3.2400 
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FIG. 22: (a) Plot of the 7r°7r° invariant mass versus the 
K + K~ invariant mass for all selected K + K~ events; (b) 

the 7r°7r° invariant mass projections for events in the (j> peak 
(open histogram), sidebands (cross-hatched), and control re- 
gion (hatched); (c) the K + K~ mass projection for events in 
the signal (open) and control (hatched) regions; (d) the dif- 
ference between the open histogram and sum of the other 
contributions to (b). 



We obtain the number of e + e~ — > </>7r°7r° events in 
0.04 GeV/c 2 intervals of 07r°7r° invariant mass by fitting 
the K + K~ invariant mass projection in that interval to 
the <j> signal, after subtracting the non-K + K~tt°tt° back- 
ground, the same way as described in Sec. IIVFI The ob- 
tained cross section is shown in Fig. [53] and is very similar 
to that obtained from the K + K~tt + tt~ final state shown 
in Fig. 1131 The errors shown reflect not only that there 
are six times fewer events, but also a much larger back- 
ground level. 



As before, we defer discussion to Sees. I VIII and IVIII1 
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FIG. 23: Cross section for the reaction e+e"-> 0(lO2O)7r°7r° 
as a function of e + e~ cm. energy obtained from the 
K+K-n°n final state. 



F. The </>(1020)/ (980) Intermediate State 

Since the background under the /o(980) peak in 
Figs. I2"27 b.d) is 25% or less, we are able to extract the 
</>(1020)/ (980) contribution. As in Sec. lTVGl we require 
the dipion mass to be in the range 0.85-1.10 GeV/c 2 and 
fit the background-subtracted K + K~ mass projection in 
each 0.04 GeV/c 2 interval of K + K~ir°ir° mass to obtain 
the number of <f>fo events. Again, some 4>tt°tt events are 
present in which the 7r°7r pair is not produced through 
the f - 

We convert the number of /o(980) events in each 
mass interval into a measurement of the e + e~ — > 
</>(1020)/o(980) cross section as described previously, and 
divide by the /o(980)— > 7r°7r° branching fraction of 1/3 
to obtain the /o(980) — > tttt value. The cross section, 
corrected for the 0(1020) — > K + K~ decay rate, is shown 
in Fig. [M] as a function of E c m . and is listed in Table llXl 
Due to the smaller number of events, we have used larger 
intervals at higher energies. The overall shape is consis- 
tent with that obtained from the K + K~tt + tt~ final state 
(see Fig. H31) . and there seems to be a sharp drop near 
2.2 GeV; however, the statistical errors are large and no 
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FIG. 24: Cross section for the reaction e + e~ — > 

0(lO2O)/o(98O), /o — > tttv as a function of e + e~ cm. energy 
obtained from the K + K~-k°it final state. 



conclusion can be drawn from this mode alone. Possible 
interpretations are discussed in Section [Villi 



Signal 
region 




Control 
region 



-i I I 



lilts 



10 



20 



30 



40 



50 60 
5C 4K 



FIG. 25: Distribution of \ 2 from the three-constraint fit 
for K + K~ K + K~ candidates in the data (points). The open 
histogram is the distribution for simulated signal events, nor- 
malized as described in the text. The shaded histogram rep- 
resents the background from non-ISR events, estimated as 
described in the text. The region defined by the dashed con- 
tour is for simulated ISR K + K~ir + ir~ events with at least 
one pion misidentified as a kaon. 



VI. THE K+R-R+R- FINAL STATE 
A. Final Selection and Background 

Figure US] shows the distribution of x\k f° r the 
K + K~K + K~ candidates as points. The open histogram 
is the distribution for simulated K + K~ K + K~ events, 
normalized to the data in the region x\k < 5 where the 
relative contributions of the backgrounds and radiative 
corrections are small. The shaded histogram represents 
the background from non-ISR e + e~ — > qq events, eval- 
uated as for the other modes. The region defined by 
the dashed contour represents the background from sim- 
ulated ISR K + K~tt + tt~ events with at least one charged 
pion misidentified as a kaon. 

We define signal and control regions by \\k < 20 and 
20 < x\k < 40, respectively. The signal region contains 
4190 data and 14904 simulated events, and the control 
region contains 877 data and 1437 simulated events. Fig- 
ure shows the K + K~ K + K~ invariant mass distri- 
bution from threshold up to 4.5 GeV/c 2 for events in 
the signal region as points with errors. The qq back- 
ground (shaded histogram) is small at all masses. Since 
the ISR K + K~tt + it~ background does not peak at low 
x\k values, we include it in the background evaluated 
from the control region, according to the method ex- 
plained in Sec. II V Al It dominates this background, which 
is about 20% for 2.3-2.6 GeV/c 2 and 10% or lower at all 
other mass values. The total background is shown as the 
hatched histogram in Fig. [211 



We subtract the sum of backgrounds from the number 
of selected events in each mass interval to obtain the num- 
ber of signal events. Considering the uncertainties in the 
cross sections for the background processes, the normal- 
ization of events in the control region, and the simulation 
statistics, we estimate that the systematic uncertainty on 
the signal yield is less than 5% in the 2-3 GeV/ c 2 region, 
but it increases to about 10% above 3 GeV/c 2 . 

B. Selection Efficiency 

The detection efficiency is determined as for the other 
two final states. Figure l2"7Ta) shows the simulated 
K + K~ K + K~ invariant- mass distributions in the signal 
and control regions from the phase space model. We di- 
vide the number of reconstructed events in each mass 
interval by the number generated in that interval to ob- 
tain the efficiency shown by the points in Fig. [27T b) . It is 
quite uniform, and we fit the measurements using a third- 
order polynomial, which we then use to obtain the cross 
section. As discussed previously, this efficiency includes 
the difference between the EMC acceptance and the re- 
gion of ISR photon simulation. A simulation assuming 
dominance of the 4>K + K~ channel, with the K + K~ pair- 
in an angular-momentum S-wave state, gives consistent 
results, as shown by the dashed curve in Fig. 127T b). and 
we estimate a 5% systematic uncertainty associated with 
the difference. We correct only for mis-modeling of the 
track-finding and ISR-photon-detection efficiency as in 
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C. Cross Section for e + e 
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FIG. 26: Invariant mass distribution for K + K~ K + K~ can- 
didates in the data (points). The shaded histogram represents 
the non-ISR background, and the hatched region is for the ISR 
background from the control region, which is dominated by 
the contribution from misidentified ISR K + K~~ -k+it~ events. 
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FIG. 27: (a) Invariant mass distributions for simu- 
lated K+ K~~ K+ K~ events in the signal (open) and control 
(hatched) regions (see Fig. I25[) : (b) net reconstruction and 
selection efficiency as a function of mass obtained from this 
simulation; the curves represent third-order polynomial fits 
for the phase space model (solid) and the <f>K+ K~ model 
(dashed). 



We calculate the e + e~— !• K+K ~K+K~ cross section 
in 0.025 GeV intervals of E c . m . from the analog of Eq.(|2]), 
using the invariant mass of the K+K~ K+K~ system to 
determine the cm. energy. We show the cross section in 
Fig. [2H1 an d list the measured values in Table [X] The 
cross section increases from threshold to a peak value of 
about 0.1 nb near 2.7 GeV, then decreases slowly with 
increasing energy. The only statistically significant nar- 
row structures are the large Jji\) peak and a possible 
narrow structure near 2.3 GeV, which will be discussed 
in Sec. IVIDI Again, the differential luminosity contribu- 
tion in each E c m . interval includes corrections for vacuum 
polarization that should be omitted for the calculations 
of a M . This measurement supersedes our previous re- 
sult ni. 




(GeV) 



Sec. UVBl 



cm. 

FIG. 28: The e + e"-s> K+K~ K+ K~ cross section as a func- 
tion of e+e~ cm. energy measured with ISR data at BABAR. 
The errors are statistical only. 



The simulated K+K~K+K~ invariant mass resolu- 
tion is 3.0 MeV/c 2 in the 2.0-2.5 GeV/c 2 range, increas- 
ing with mass to 4.7 MeV/c 2 in the 2.5-3.5 GeV/c 2 range, 
and to about 6.5 MeV/c 2 at higher masses. Since the 
cross section has no sharp structure except for the J /%p 
peak, we again make no correction for resolution. The 
errors shown in Fig. [2H an d listed in Table |X] are sta- 
tistical only. The point-to-point systematic uncertainties 
are much smaller, and the errors on the normalization 
are summarized in Table IXI1 along with the corrections 
applied to the measurements. The total correction is 
+4.0%, and the total systematic uncertainty is 9% at 
low mass, linearly increasing to 13% above 3 GeV/c 2 . 
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TABLE X: Summary of the cross section measurements for e + e — > K + K K + K . Errors are statistical only. 



E c . m . (GeV) 
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(nb) 


E c . m 


. (GeV) 
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(nb) 


Eo.m. (GeV) 
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(n 


b) 


Earn. (GeV) 


a 


(n 


b) 


2.0125 


0.002 


± 


0.002 


2 


6375 


0.100 


± 


0.016 


3.2625 


0.035 


± 


0.008 


3.8875 


0.020 


± 


0.005 


2.0375 


0.003 


± 


0.004 


2 


6625 


0.083 


± 


0.013 


3.2875 


0.030 


± 


0.009 


3.9125 


0.011 


± 


0.005 


2.0625 
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FIG. 29: Invariant mass distribution for all K + K~ pairs in 
selected e + e~ K + K~ K + K~ events (open histogram), and 
for the combination in each event closest to the 0-meson mass 
(hatched). 



D. The 4>{1020)K + K- Intermediate State 

Figure [2U shows the invariant mass distribution for 
all K + K~ pairs in the selected K + K~ K + K~ events 
(4 entries per event) as the open histogram. A promi- 



TABLE XI: Summary of corrections and systematic uncer- 
tainties for the e + e~ —¥ K + K~ K + K~ cross section measure- 
ments. The total correction is the linear sum of the individual 
corrections, and the total uncertainty is the sum in quadra- 
ture of the separate uncertainties. 



Source 


Correction 


Uncertainty 




Rad. Corrections 




1% 






Backgrounds 




5% , E c . m . 


< 3 


GeV 






5-10% , E c . m . 


> 3 


GeV 


Model Dependence 




5% 






X4K Distribution 




3% 






Tracking Efficiency 


+3.0% 


2% 






Kaon ID Efficiency 




4% 






ISR-photon Efficiency 


+ 1.0% 


0.5% 






ISR Luminosity 




3% 






Total 


+4.0% 


9% , E c . m . 


< 3 


GeV 






9-13% , Ec.m. 


> 3 


GeV 



nent (f> peak is visible along with a possible excess near 
1.5 GeV/c 2 . The hatched histogram is for the pair in 
each event with mass closest to the nominal <j) mass, 
and indicates that the <pK + K~ channel dominates the 
K + K~ K + K~ final state; we do not see any other sig- 
nificant contribution. If the invariant mass of the K + K~ 
pair that is closest to the <f> mass is within ±10 MeV/c 2 
of the <f> peak, then we include the invariant mass of the 
other K + K~ combination in Fig. l5OT a). Events with 
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FIG. 30: (a) The invariant mass distribution for K + K~ pairs in events in which the other K + K~ pair has mass closest to, and 
within 10 MeV/c 2 of, the nominal <f> mass (open histogram); events within ±50MeV/c 2 of the J/ip mass have been excluded. 
The hatched histogram corresponds to events with K + K~ K + K~ invariant mass in the J/tp peak. The numbered regions of 
the combined histograms from (a) are used to calculate the cross sections shown in Figs, (b), (c) and (d) for regions 1, 2 and 
3 respectively. 



K + K~ K + K~ mass within ±50 MeV/c 2 of the J/i/j mass 
are excluded. Events within ±50 MeV/c 2 of the J/tp 
mass are shown as the hatched histogram. The latter is 
in agreement with results from the BES experiment [3l[ , 
for which the structures around 1.5, 1.7, and 2.0 GeV/c 2 
were studied in detail. For the dots with error bars there 
is an enhancement at threshold that can be interpreted 
as being due to /o(980) —> K + K~ decay. This is ex- 
pected in light of the <fifo cross sections measured above 
in the K + K~tt + tt" and K + K~ir°ir° final states, but a 
contribution from the ao(980) — > K + K~ can not be ex- 
cluded. For the combined histograms of Fig. ISPTa). we 
select events with m(K + K~) < 1.06 GeV/c 2 (shown as 
region 1) and calculate a cross section enriched in the 
e+e- ->■ 0/ o (98O) reaction (Fig. GSJb)). A bump at 
E c .m. = 2.175 GeV is seen; however, the small number 
of events and uncertainties in the /o(980) — > K + K~ line- 
shape do not allow a meaningful extraction of the cross 
section for this /o(980) decay mode. 

A clear signal corresponding to the /2(1525) is seen 
in both histograms shown in Fig. 1307 a). The /2(1525) 
region is defined by 1.45 < m(K + K~) < 1.6 GeV/c 2 , and 
is indicated as region 3 in Fig. [307a). The corresponding 
cross section is shown in Fig. l30T d) and exhibits a broad 
(about 0.10-0.15 GeV) structure at 2.7 GeV and a strong 
J/ijj signal. In Fig. I3"07a)fopen histogram) there is an 
indication of structure for K + K~ invariant mass in the 
1.3-1.4 GeV/c 2 region; this may correspond to production 
of the 0/ o (137O) final state. 

Finally, we tried to find a region of K + K~ invariant 
mass corresponding to the spike seen at about 2.3 GeV 
in the total e + e~ — > K + K~ K + K~ cross section shown 
in Fig. 1281 This spike is much more significant if we 
require 1.06 < m(K + K~) < 1.2 GeV/c 2 , shown as region 
2 in Fig. IBUT a) , with corresponding cross section shown 
in Fig. l30T c). We have no explanation of this structure. 

We observe no significant structure in the K + K~ if ± 
mass distribution. 

We use the cf>K + K~ events to investigate the possi- 
bility that part of our <fiir + ir- signal is due to cj>K + K~ 



events with the two kaons interpreted as pions. No struc- 
ture is present in the resulting K + K~tt + it~ invariant 
mass distribution. 



VII. THE e+e~ </>7T7r CROSS SECTION 
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FIG. 31: The two-Breit-Wigner fit to the n tt~ invariant 
mass distribution of Fig. ll2[ d). The dashed curve corresponds 
to the inclusion of the partial width to KK in the propagator 
of the / (980) BW. 

We next perform a more detailed study in the E c m . re- 
gion from threshold to 3.0 GeV of the e+e" -> 0(1O2O)7T7t 
cross section. For this study we use the cross section for 
the </>7r + 7r - final state shown in Figs. [TBI and [T5l after scal- 
ing by a factor of 1.5 to take into account the 4>tt°tt° con- 
tribution. The cross section for the 4>ir°ir° (see Fig. |2"3")) 
final state does not help much due to large statistical er- 
rors. There are at least two candidate resonant structures 
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in Fig. [13] These are associated with the peaks observed 
at 1.7 GeV and at 2.1 GeV. As shown in Sec HYG] the 
latter is related to </>(1020)/o(980) production, while the 
best candidate for the former may be the </>(1680), which 
is a radial excitation of the ss state decaying predomi- 
nantly to K*(892)K [ll]. This would be another confir- 
mation of the decay of this state to </>(1020)7r7r, previously 
reported in Refs. 0,0. 

As discussed in Sec. lIVFl we associate the narrow peak 
in the 7r + 7r~ invariant mass distribution, shown on a 
larger scale in Fig. [3TJ with the / (980) (denoted as the 
fo meson), and observe a broad enhancement at about 
0.6 GeV/c 2 ; the angular distributions of Fig. [14] justify 
that these structures are in an S-wave state. This low 
mass bump can not be formed by pure three-body phase 
space. Indeed, the (/>(1020)7T7r threshold is 1.3 GeV, 
but the observed cross section has a slow rise starting 
at 1.4 GeV. This indicates that the observed structure 
could be a result of /o(600) resonance decay. 

The observed two-pion-mass shape of the /o(600) (de- 
noted as the a meson) is distorted by the 0(1O2O)7T7t final 
state. This is less of an issue for the narrower /o(980). 
Nevertheless, to obtain mass and width parameter values 
for these states, we fit the data of Fig. [31] using a function 
consisting of an incoherent sum of two S-wave relativistic 
BW intensity distributions, modified to account for the 
two pion phase space. The fit values obtained are 

m„ = (0.692 ±0.030) GeV/c 2 ,^ = (0.538 ±0.075) GeV, 

(3) 

and 

m fo = (0.972±0.002) GeV/c 2 ,r /o = (0.056±0.011) GeV, 

(4) 

and the fit result is represented by the solid curve in 
Fig. [SB Note that the /o(980) parameters are consistent 
with the PDG values [5] , indicating that interference with 
the /o (600) (or tttt coherent continuum) is minimal. This 
is expected because events with m(Tnr) < 0.85 GeV/c 2 
are associated with the resonance at 1.7 GeV/c 2 in the 
</>(1020)7T7r mass, while the /o(980) contributes only to 
a structure above 2 GeV/c 2 (see Fig. [T5]) . To confirm 
this, we examine two m(7r + 7r~) distributions using the 
selections shown in Fig. [12] but for events with either 
m(07r+7r~) < 1.95 GeV/c 2 or 1.95 < m(07r + 7r~) < 
3.0 GeV/c 2 . For the first case we observe only the bump 
at 0.6 GeV/c 2 of Fig. [3TJ with no evidence for the / (980). 
For the second case we see a clear / (980) signal but no 
evidence for the / (600). We fit each distribution the 
same way as the data in Fig. 1311 The resulting param- 
eters for the /o(600) and /o(980) are in agreement with 
those presented above. 

The dashed curve of Fig. (3JJ is obtained when the 
/o(980) — > KK partial width is incorporated into the 
BW propagator (the so called Flatte approximation used 
in Ref. [301 ] with parameters cl/c2 and m ■ cl, which cor- 
respond to the ratio of the coupling constants g 2 KK lg\ v 
and effective /o(980) width). It differs only slightly at 
the top of the /o(980), but the wider shape of the Flatte 



function leaves less room for the remaining events and 
we obtain: 

m a = (0.631 ±0.020) GeV/c 2 ,r CT = (0.472 ±0.075) GeV. 

The obtained Flatte function parameters are in agree- 
ment with those obtained in Ref. [3(|: c2/cl = 2.20 ± 
0.67, in • cl — 0.131 ± 0.033. 

The Flatte approximation gives a little better descrip- 
tion of the observed tttt mass spectrum, and so we use 
it in the analysis of the structures observed in the 4>irir 
cross section. 

It appears that the structure at E c m . m 2.1 GeV in 
the (f>ir + TT~ cross section fFigs. [TBI and IT51) couples to the 
/o(980) but not to the Jo (600). This is very similar to the 
behavior observed for the tt + tt~ system in J/ip — > </)7r + 7r~ 
decay [3l[ (and demonstrated with our data in Fig. |4"21 of 
Sec. Iixp . and in D s —t tt + tt~tt + decay (32j. In both in- 
stances a clear Jo (980) signal is observed, while the broad 
/o(600) enhancement of Fig.[3TJis absent. In contrast we 
note that in J/tp — > ujtt + tt~ decay [33| exactly the op- 
posite behavior is observed; the tt + tt~ system exhibits a 
broad low-mass enhancement, and there is no evidence 
of an /o(980) signal. 

In contrast with the "clean" m{<j)TT + TT~) distribution, 
obtained from the fit on the <f> peak, the m(rT + TT~) dis- 
tribution is obtained by the selection of <j> signal in the 
K + K~ invariant mass distribution, with background 
subtraction performed using the (f> side bands and control 
region of the x 2 distribution (see Fig. IT2"]) . To minimize 
these uncertainties, we use a BW description for a and 
the Flatte approximation for / to incorporate these two 
states in a simple model describing the structures in the 
(fmw cross section data of Fig. [T3]( after scaling by a factor 
of 1.5 to take into account the </>7r°7r° contribution). The 
model consists of the incoherent addition of two contri- 
butions at each value of E c m .. The first represents the 
decay process </>(1680) —> (/>/o(600), with the parame- 
ters of the a given by Eq. ([5]); the second results from 
the coherent superposition of amplitudes describing the 
processes </>(1680) <£/ (980) and V(2175) -> </>/ (980), 
where the Y(2175) BW amplitude describes the peak ob- 
served at w2.2 GeV in Fig.[U We note that in Ref. [|] 
the contribution from </>(1680) — > </>/o(980) decay was not 
taken into account. We see no physical evidence to justify 
doing this, and so allow the presence of this amplitude 
in our model. The angular distributions of Fig. [14] are 
consistent with the </>(1020) and the S-wave tttt system 
being in an S-wave orbital angular momentum state, and 
so our model includes no centrifugal barrier factor in the 
amplitude representations. 

We fit the observed e + e~ — > (f)inr cross section using 
the function 



iV(s) 
s 3 / 2 

P<t>fa ( s ) 
s 3/2 



Au(s) 



(6) 



A 22 {s) 
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FIG. 32: (a) The fit to the e + e~ —¥ 07T7T cross section using the model described in the text; the entire contribution due to the 
0(1680) is shown by the dashed curve. The dotted curve shows the contribution for only 0/o(98O) decay, (b) Comparison of 
the data and the curve obtained from the overall fit, with the restriction m(7T7r) < 0.85 GeV/c 2 . (c) The e + e~— >0(lO2O)/o(98O) 
cross section with the requirement 0.85 < m(Tnv) < 1.1 GeV/c 2 ; the dashed and dotted curves represent the contributions from 



0(1680) -> 0( 1020) /o (980), and 0(1680) 
cross section data. 



0( 1020) /o (600) calculated using the parameter values from the overall fit to the 



TABLE XII: Summary of parameter values obtained from the fits with Eq. ([SJl described in the text. An asterisk denotes a 
value that was fixed in that fit. 



Fit 


All rn(mv) 


m(7T7r) < 0.85 GeV/c 2 


0.85 < m(TTiv) < 1.1 GeV/c 2 


(Til (nb) 
mi(GeV/c 2 ) 


0.655±0.039±0.040 


0.678±0.047±0.040 


0.655* 


1.742±0.013±0.012 


1.733±0.010±0.010 


1.742* 


ri(GeV) 


0.337±0.043±0.061 


0.300±0.015±0.037 


0.337* 


(J22 (nb) 
m 2 (GeV/c 2 ) 


0.082±0.024±0.010 


0.082* 


0.094±0.023±0.010 


2.176±0.014±0.004 


2.176* 


2.172±0.010±0.008 


r 2 (GeV) 


0.090±0.022±0.010 


0.090* 


0.096±0.019±0.012 


cri2(nb) 


0.152±0.034±0.040 


0.152* 


0.132±0.010±0.010 


tp (rad) 


-1.94±0.34±0.10 


-1.94* 


-1.92±0.24±0.12 


X 2 /n.d.f. 


48/(67-9) 


46/(66-4) 


38/(46-6) 


P(x 2 ) 


0.74 


0.96 


0.40 



variant mass available at E c m . = y/s, and are obtained 
from 

P$nn{s) = / BWT V7T (m)q(s,in,m^)dm, (7) 

J2m 7r 



where 



with 



, 3/2 -p 



mf ~ s ~ iy/sTi(s) ' 



i = 1 for the 0(1680), i = 2 for the Y(2175), 
j = 1 for the / (600), j = 2 for the /q(980), 

so that 

An(s) describes 0(1680) 0(1O2O)/ O (6OO) decay, 
A 12 (s) describes 0(1680) 0(1O2O)/ O (98O) decay, 
A 22 (s) describes F(2175) -> 0(1O2O)/ O (98O) decay; 

s = m , mi and Ti are the mass and width of 
the 0(1680), 77i2 and are the mass and width of the 
y(2175), and the cry represent the peak cross section 
values. 

The factors P < p (y {s) and P<pf (s) represent quasi- two 
body phase space integrated over the range of 7T7t in- 



where BW ww (m) is a BW function with /o(980) param- 
eters (BWf (m)) to define P<j>f (s), or with /o(600) pa- 
rameters (BW a (m)) to define P<p a (s) and q is the 
momentum of the particles with masses m and in the 
two-body reaction at E c . m . = y/s. 

Since the decay 0(1680) -> 0(1O2O)/ O (98O) is sup- 
pressed by phase space near ^/s = mi, the value of a\ 2 
is much smaller than that of cth, but its contribution to 
a(s) increases rapidly beyond the 0(1O2O)/ O (98O) thresh- 
old. 

The 0(1680) resonance decays mainly to KK* (892) 
and 0(1020)7/ 0, B3]. We find that it has a branching 
fraction of about 10% to 07T7t, which together with other 
modes listed in PDG, leads to an energy-dependent width 
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that can be written as 



ri(s) 



0.2 



0.7 



m\P 2K (s) 



0.1 



. (8) 



with P 2 k(s) = t? \\fs-, m K ,m K *), and i%(s) = 
q(y/~8,m<f l ,m 7} ). 

For the second resonance candidate, which decays 
mostly to </f>7r7r, the energy dependence of the width is 
written as 



r a (*) 



m 2 P <t> ^(s) 
3 V2p 07r7r ( m 2) 



(9) 



We note that the introduction of an energy dependence 
for each width significantly increases the values of the res- 
onance mass and width, especially for broad structures. 

The results of the fits are shown in Fig. [32] an d sum- 
marized in Table IXII1 The first error is statistical, and 
the second error represents the systematic uncertainty 
estimated as a difference in fitted values for two differ- 
ent descriptions of the two-pion spectrum as shown in 
Fig. [3T] In Fig. 1327a) we show the contribution from 
the (£(1680) for both modes (dashed curves), and for 
(£(1680) <£(1020)/o(980) only (dotted curve). The in- 
crease of the cross section at about 2 GeV is explained by 
the opening of the (£/o(980) decay channel of the (£(1680) 
resonance. However the fit shows that an additional rela- 
tively narrow state is needed in order to provide a better 
description of the observed data. 

It is important to note that this model describes the 
observed data very well independently of the m(7T7r) re- 
gion selected. Figure l32T b) shows the ifiirw cross section 
for m(7r + 7r~) < 0.85 GeV/c 2 for the data; the curve is 
obtained by using the parameter values from the overall 
fit and yields x7 n -d.f. = 63/(66-1) (P(x 2 ) = 0.54). If 
we fit this distribution, slightly better parameter values 
can be obtained (see Table IXlT)) . but these still agree well 
with those from the overall fit. We consider them as our 
measurement of the (£(1680) resonance parameters. They 
correspond to the product of the electronic width, T ee , 
and branching fraction to 4>tttt, B^tt, 



. r — 

P.P. 



12ttC 



(42 ± 2 ± 3) eV , 



where we fit the product TiCTh to reduce correlations, 
and C, the conversion constant, is 0.389 mb( GeV/ c 2 ) 2 [f| . 
The second error is systematic, and corresponds to the 
normalization uncertainty on the cross section, and to 
the uncertainty in the m(7T7r) distribution description. 

If we require 0.85 < m(nn) < 1.1 GeV/c 2 (Fig. E^c)), 
then without additional fitting the model yields x 2 / n -d.f. 
= 48/(46-1) (P(x 2 ) = 0.31), and improves to x 2 /n.d.f. = 
38/(46-6) (P(x 2 ) = 0.40 ) by refitting using the param- 
eter values listed in Table IXII1 If we try to explain the 
observed cross section only in terms of the (£(1680) with- 
out any narrow state (dashed curve in Fig.[3Stc)), the fit 



gives x 2 /n.d.f. = 123/(46-2) (P(x 2 ) = 10~ 7 ) and so this 
hypothesis is not compatible with the data. Note, that 
the contribution of (£/o(600), shown by dotted curve in 
Fig-EUJc), is very small. 

The model described above provides an excellent de- 
scription of the observed cross section behavior, and sug- 
gests that the V(2175) may not be a radially excited 
ss state, since such a state would be expected to be 
much wider (300-400 GeV/c 2 ) and also should decay to 
(£/o(600), like the (£(1680). 



VIII. 



<£/o NEAR THRESHOLD 



The behavior of the e + e~ — > (£/o cross section near 
threshold shows a structure near 2.175 GeV, and we have 
published this result in Ref. 7]. Here we provide a more 
detailed study of the cross section for this channel in 
the 1.8-3 GeV region with the full BABAR dataset. In 
Fig. [33J we superimpose the cross sections measured in 
the K + K~ii + ir~ and K + K~ir a ir® final states (shown in 



Figs . [151 and [24]) : they are consistent with each other. 

We perform a combined fit to these cross section data 
using Eq.© with the two-pion mass restricted to the re- 
gion 0.85-1.1 GeV/c 2 . We fix the (£(1680) parameters 
for the (£(1020)/o(600) decay mode (which gives a small 
contribution in this mass range) and allow all other pa- 
rameters to float. The result of the fit is shown as the 
solid curve in Fig. [33] As demonstrated in Ref. Q, 
the observed pattern can be a result of a constructive 
or destructive interference of the narrow structure at 
2.175 GeV with the coherent background. 

The fit with constructive interference gives the reso- 
nance parameter values 

cr 22 = (0.093 ± 0.021 ± 0.010) nb, 

m 2 = (2.180 ±0.008 ±0.008) GeV/c 2 , 

T 2 = (0.077 ±0.015 ±0.010) GeV, 

<£ 2 = (-2.11 ± 0.24 ± 0.12) rad, 

cr 12 = (0.140 ± 0.009 ± 0.010) nb, 

andx 2 /n.d.f.= 57/(61-6) (P(x 2 ) = 0.33). The statistical 
precision is improved compared to that of Ref. 0], for 
which the analysis was based on half as much data. For 
this state we estimate the product of electronic width 
and branching fraction to (£/o as 



B 



4>fo 



r — 
x pp. — 



r 2 er 22 TO 2 

12ttC 



(2.3 ±0.3 ±0.3) eV 



where we fit the product r 2 <r 22 to reduce correlations. 
The second error is systematic, and corresponds to the 
normalization uncertainty on the cross section. 

The destructive interference yields exactly the same 
overall curve with the same parameters for the mass 
and width of the narrow state, but significantly larger 
peak cross section with opposite sign of the mixing angle: 
cr 22 = (1.13±0.15±0.12)nb, ip 2 = (2.47±0.17±0.13)rad. 
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FIG. 33: The e+e" 4>{ 1020) f (980) cross section mea- 
sured in the K + K~ tt + it~ (solid dots) and K + K~ ir°ir° (open 
squares) final states. The solid (dashed) curve represents 
the result of the two-resonance (one-resonance - 0(1680) — ► 
4>{ 1020) f (980)) fit using Eq.@, as described in the text. The 
hatched area and dotted curve show the Y(2175) contribution 
for two solutions. 




E c . m . (GeV) 

FIG. 34: The e+e"->0(lO2O)/o(98O) cross section measure- 
ments from the K + K~ tv + -k~ final state from BABAR (dots) 
and Belle [j| (squares) . 



shown as the dashed curve in Fig. [33] It is a poor fit to 
the region below 2.3 GeV, but gives a good description of 
the cross section behavior at higher values of E c m .. The 
fit, with or without the resonance at 2.18 GeV/c 2 , gives a 
maximum value of the 0(1680) — > 0/o cross section of 0.3 
nb at E c m . w 2.1 GeV. This is of independent theoretical 
interest, because it can be related to the </> — > /o(980)7 
decay studied at the 0- factory [28|, |29j . 

The significance of the structure calculated from the 
change in \ 2 between the fits with and without the res- 
onance at 2.18 GeV is V150 — 61 = 9.4 standard devi- 
ations; the x 2 value, 61 for 61-2 n.d.f., yields the same 
probability as the \ 2 value 57 for 61-6 n.d.f.. 

The cross section measurements from the K + K~n + ir~ 
final state shown in Fig. EEU are compared to those from 
Belle [9( in Fig. [34] There is good overall agreement 
between the results from the two experiments. Over- 
all agreement between the results of the fits to the 
BABAR and Belle data is also good. 

A. Structures in the K + K~f (980) final state 

We next search for other decay modes of the V(2175) 
state. Figure l3"5T a) shows the "raw" (no background 
subtraction) two-pion mass distribution for all selected 
K + K~ir + Tr~ events, and Fig. IBl^T a) shows the same dis- 
tribution for the K+K-TT 7r° sample. The / (980) con- 
tribution is relatively small for the charged-pion mode, 
and larger for the neutral-pion mode. If we select the 
region 0.85 < m(Tnr) < 1.1 GeV/c 2 , and plot the 
K + K~irir mass distribution, the bump at 2.175 GeV/c 2 
is seen much more clearly in spite of larger background 
(Figs. ESJb) and |3gfb)), and a bump at 2.5 GeV/c 2 
is also seen; the rest of events have no structures at 
2.175 GeV/c 2 or 2.5 GeV/c 2 (Figs. E^b) and H^b) 
hatched histograms). The bumps are seen only in the 
K+K- / (980) sample (Figs. |5STc) and He)), but if 
we select the 0(1020) region, no bumps are seen at 
2.5 GeV/c 2 , as shown by the hatched histograms in 
Figs. [3itc) and EStc). 

From these histograms we can conclude that the 
V(2175) resonance has a K + K~ /q(980) decay mode 
when the K + K~ system is not from 0, and that the 
decay rate is comparable to that for 0/o(98O). Also an- 
other state at 2.5 GeV seems to exist; this decays to 
K + K~ / (980) (but seems not to couple to 0/ o (98O)) 
with width w0.06-0.08 GeV (see Ref. @|). The large 
background does not allow us to clearly separate this 
state. 



To select between two solutions, we need more informa- 
tion on the decay rates to another modes, which are not 
available now. 

If we assume no resonance structure other than the tail 
from 0(1680) ->■ 0(1O2O)/ O (98O), the fit yields x 2 /n.d.f.= 
150/(61 - 2) with P(x 2 )=8 ■ 10" 9 . The result of this fit is 



IX. THE CHARMONIUM REGION 

For the E c . m . region above 3 GeV, our data can be used 
to measure, or set limits on, the decay branching frac- 
tions for the J/tp and ^(2S) (See Figs. H \M and [25J). 
In addition, these signals allow checks of our mass scale, 
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FIG. 35: (a) The m(7r + 7r~) distribution without background subtraction for K + K~ tt + tt~ events. The vertical lines indicate 
the /o(980) region, (b) All selected K + K~ tt + -k~ events (open histogram), selected K + K~ fo(980) events (cross-hatched 
histogram), and all the rest (hatched histogram), (c) The K + K~ /o(980) events (open histogram) in comparison with the 
4>{ 1020) fo (980) sample (hatched histogram). 
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FIG. 36: (a) The m(7r°7r°) distribution without background subtraction for K + K~ -7r°7r° events. The vertical lines indicate the 
/o(980) region, (b) All selected K + K~ tt°iy events (open histogram), selected K + K~ fo(98Q) events (cross-hatched histogram), 
and all the rest (hatched histogram), (c) The K + K~ fo(980) events (open histogram) in comparison with the 0(1020) fo (980) 
sample (hatched histogram). 



and of our measurements of mass resolution. Figure [37] 
shows the invariant mass distributions for the selected 
K+K-tt+tt-, K + K-T7 Q n°, and K + K~K + K~ events, 
respectively, in this region, using smaller mass intervals 
than in the corresponding Figs. [2] [TTj and [26] We do 
not subtract any background from the K + K~ir + ir~ and 
K + K ~ K + K~ distributions, since it is small and nearly 
uniformly distributed, but we use the x\k2-k° con trol re- 
gion to subtract part of the ISR background from the 
K + K~Tr°ir° distribution. Production of the J/ip is ap- 
parent in all three distributions, and a small, but clear, 
ip(2S) signal is visible in the K + K ~-k + tt~ mode. 

We fit each of these distributions using a sum of two 
Gaussian functions to describe the J/ip signal and in- 
corporate a similar representation of a ip(2S) signal, 
although there is no clear evidence of the latter in 
Figs. 157Tb) and |57Tc). In each second-order- 
polynomial function is used to describe the remainder 
of the distribution. We take the signal function param- 



eter values from simulation, but let the overall mean 
and width values vary in the fits, together with the co- 
efficients of the polynomial. For the K + K~ir°ir° and 
K + K~ K + K~ modes we fix the ip(2S) mass position 
and take the width from MC simulation. The fits arc 
of good quality, and are shown by the curves in Fig. [57] 
In all cases, the fitted mean value is within 1 MeV/c 2 
of the nominal J/tp or iJj(2S) mass position [5( and the 
width is within 10% of the simulated resolution discussed 
in Sees. HVCl NTH and rVTCl 

The fitted J/ip signals for the K + K~tt + it~ , 
K+R-^n , and K + K~K + K~ final states are found to 
contain 3137 ± 67, 388 ± 28, and 287 ± 24 events, respec- 
tively. From the number of events in each final state /, 
Nj/^f, we calculate the product of the J/ip branching 
fraction to / and the J/tp electronic width using 



W ee 6^ 2 ■ dC/dE ■ e f (m JN ,) ■ C 



(10) 
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FIG. 37: 

e + e~ — > 
text. 



Raw invariant mass distribution for all selected events in the charmonium region for (a) e + e —¥ K + K tt + -k , (b) 
and (c) e + e~ — > K + K~ K + K~ ; in each figure the curve represents the result of the fit described in the 




1.5 1.75 2 



3.5 

2, 



2.25 2.5 2.75 3 3.25 
m(K + K"7t + 7tl(GeV/c") 

FIG. 38: The K it^ invariant mass versus K + K~tt + h~ 
invariant mass for events with the K^n^ combination in one 
of the 7^*(892)° regions of Fig. [5ja); for events in overlap 
region, only one combination is chosen. 
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FIG. 39: The K ± ty^ 1 mass projection for events from Fig. 1381 
with K + K~ -k + it~ invariant mass within 50 MeV/c 2 of the 
nominal J/ip mass (open histogram), and for events for which 
this mass value is 50-100 MeV/c 2 less than nominal (hatched). 



where dC/dE = 173.1 ± 1.7 nb~7MeV, and e f {mj/^) 
are the ISR luminosity and corrected selection efficiency, 
respectively, at the J/tp mass, and C is the conver- 
sion constant. We estimate £k+k--k+-k- = 0.198 ± 
0.006, e K + K - w o w o = 0.079 ± 0.004, and e K+K - K+K - = 
0.173 ± 0.012 using the corrections and errors discussed 
in Sees. HVCl IVCl and lyiCl 

We list the values of the product of the branching 
fraction(s) and T J e i^ in Table IXIII1 and using Y J J e ^ — 
(5.55 ± 0.14) keV [5j, obtain the corresponding branch- 
ing fraction values and list them together with their PDG 
values [f| ■ The systematic uncertainties quoted include a 

2.5% uncertainty on vii^ . Our measured branching frac- 
tions of K+K-it+TT', K+K-n^n", and K + K~ K + K~ 



are more precise than the current PDG values, which 
are dominated by our previous results ((6.6±0.5) x 10~ 3 , 
(2.5±0.3)xlCT 3 and (7.6±0.9) x 1(T 4 , respectively Q). 

These fits also yield 133±21 K+R-tt+tt- events, 17±9 
K + K-T7°Tr events and 13±6 K + K~K + K~ events in 
the ■0(25') peak. We expect 12 events from •0(25') — > 
J/ipn + Tr~ — > K + K~tt + tt~ from the relevant branching 
fractions ||, which is less than the statistical error. Sub- 
tracting this contribution and using the calculation anal- 
ogous to Eq.|l0]), with dC/dE = 221.2 ±2.2 nb^/MeV, 
we obtain the product of the branching fraction and 
electronic width for the decays ifj^S) — > K + K~ir + ir~ , 
tp(2S) ->■ K + K-n°n°, and ^{2S) -> K+R-K+R-. 
Dividing by rt (2S) = 2.36 ± 0.04 keV @, we ob- 
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FIG. 40: Raw invariant mass distributions in the charmo- 
nium region for (a) candidate e + e~— > 4>it + -k~ events (open his- 
togram), and for events in the cj> sideband regions of Fig.[l2jc) 
(hatched); (b) candidate e + e~ — > 4>-k°tt events (open his- 
togram) and events in the X2K2-n° control region (hatched). 
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FIG. 41: Raw invariant mass distribution in the charmo- 
nium region (a) for candidate <pfo, fo— >-k + tt~ events (open 
histogram), and for events in the <f> sideband region (hatched), 
and (b) for candidate 4>fo, /o^7r°7r° events (open histogram) 
and for events in the x\k2-k° control region (hatched). 



tain the branching fractions listed in Table IXIIII The 
K + K~iy + tt~ and # + #~# + #~ values are consistent 
with those in Ref. joj . There is no entry in Ref. Q for 
the K+K-n^n" deca y mod e of the tp{2S). 

As noted in Sec. IIVDI and shown in Fig. [5] and 
Fig. [51 the K + K~n + ir~ final state is dominated by the 
#*(892)°#~7r + channels, with a small contribution from 
the #*(892)°#*(1430)° channels. Figure |38] shows a 
plot of the invariant mass of a K tt^ pair versus that 
of the K + K~tt + it~ system for events with the mass of 
the # T 7r ± pair near the #* (892)° mass, i.e., within the 
bands in Fig. EJa), but with only one combination plot- 
ted in the overlap region. There is a large concentration 
of entries in the J/tp band with # ± 7r =F mass values near 
1.43 GeV/c 2 , but a relatively small number of events in a 
horizontal band corresponding to the #2(1430)° produc- 
tion outside the J/tp region. We show the # ± 7r =F mass 
projection for the subset of events with #+# ~tt + it~ 
mass within 50 MeV/c 2 of the nominal J/tp mass in 
Fig. 131] as the open histogram. The hatched histogram is 
the projection for events with a K + K~ir + ir~ mass be- 
tween 50 and 100 MeV/c 2 away from the nominal J/tp 
mass. 

The Ktt distribution from the J/tp is dominated 
by the #1(1430)° and # *(1430)° signals @, HJ. A 
small signal at the #*(892)° indicates the presence of 
#*(892)°#*(892)° decay of the J/tp; this is also seen 
as an enhancement in the J/tp band in Fig. 1381 The 
enhancement at 1.9 GeV/c 2 of Fig. may be due 
to the 3 i 7 2 ground state, or to the first radial excita- 
tion of the iQ(1430), neither of which has been re- 
ported previously. Subtracting the number of side- 
band events from the number in the J/tp mass win- 
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FIG. 42: The 7r + 7r~ invariant mass distribution for 4>n + n~ 
events from the J/tp peak of Fig. 140( a) (open histogram), and 
for events in the <j> sideband region (hatched). 



dow, we obtain 710±30 events with # ± 7r =F mass in 
the range 1.2-1.7 GeV/c 2 , which we take as a measure 
of J/tp decay into #*(892)°#S, 2 (1430)°. According to 
Ref. [HI, there is an equal contribution from #q(1430)° 
and #|(1430)°, which we cannot separate with our se- 
lection. We obtain 47_± 12 events in the 0.8-1.0 GeV/c 2 
window for #*(892)°#*(892)° decay, and 185±21 events 
for decay to #*(892)°#-tt + with m(Kn) in the 1.7- 
2.0 GeV/c 2 region. We convert these to branching frac- 
tions using Eq. (jlO[) , and divide by the known branch- 
ing fractions of the #* states [5j. The results are listed 
in Table IXIIII which are more precise than those in 
Ref. Q. For the 1.7-2.0 GeV/c 2 mass region we obtain 

(2.24 ± 0.25 ± 0.15) eV. 



tU*b 



7/i/>-s--ft"*(892)°if-7r+ 



We study decays into 



and 



using the 



mass distributions shown in Figs. I4"07 a).(b). The open 
histograms are for events with K + K~ mass within the 
(p bands of Figs. H^c) and H^c). The hatched his- 
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TABLE XIII: Summary of the J/tp and tp(2S) parameters obtained in this analysis. 



ee 
ee 



v J/i> 

-L ee 

r J/v> 

-L ee 

r J/i> 

J- ee 
-L ee 
J- ee 
-L ee 
^ ee 



pV(2S) 

pV-(2S) 

r V»(2S) 
J- ee 
pV(2S) 

p-0C2S) 
1 ee 



Measured 
Quantity 



Measured 
Value ( eV) 



J/tp or i/>(2S) Branching Fraction (10 3 ) 
This work PDG2010 



™J/i)>^K+K-Tr+-n— 
"j/tp-tK+K-K+K- 



Bj/ili-Kfnr+n- ' "</)-)■ iC+Jf~ 
&J/ili-xf>K+ K~ ' &cf>^K+K 



37.94±0.81±1.10 
11.75±0.81±0.90 
4.00±0.33±0.29 



6.84±0.15±0.27 
2.12±0.15±0.18 
0.72±0.06±0.05 



•%*o 

n 0,2 



■ K-1T+ 



■B. 



■ 5 

■ B 



fx — yn~*~iv 



8.59±0. 
0.57±0. 
2.19±0. 
1.36±0. 
2.26±0. 
0.69±0. 
0.48±0. 
0.74±0. 



36±0.27 
15±0.03 
23±0.07 
27±0.07 
26±0.16 
11±0.05 
12±0.05 
12±0.05 



6.98±0. 
0.23±0. 

o.ai±o. 

0.50±0. 
1.66±0. 
0.25±0. 
0.18±0. 
0.27±0. 



29±0.21 
06±0.01 
08±0.03 
10±0.03 
19±0.12 
04±0.02 
04±0.02 
04±0.02 



V>(2S)->-K" + K ~ tt+tt- 
ip(2S)->K+K-Tr°Tr 



&i/i(2S)-*K+K- K +K- 



B 



i/>(2S)->0/ o ' B<f,-+K+K 



■ B 



fo- 



1.92±0.30±0.06 
0.60±0.31±0.03 
0.22±0.10±0.02 
0.27±0.09±0.02 
0.17±0.06±0.02 



0.81±0.13±0.03 
0.25±0.13±0.02 
0.09±0.04±0.01 
0.23±0.08±0.01 
0.15±0.05±0.01 



<j> is selected as - m(K+K~)\ < 10 MeV, B 



14, 

J/i/>-Kl>K+K- 



ob- 



taincd as 2 ■ B 

b Not corrected for the /o 
m(7r°7r°) < 1.1 GeV/c 2 

c Not corrected for the fo - 
m(ir+TT-) < 1.1 GeV/c 2 

''We compare our <j>f x , f x — > 
m(w+n-) < 1.5 with <£/ 2 (1270). 



by 0.85 < 
by 0.85 < 
mode, selected by 1.1 < 



► 7r°7r° mode, /q selected 
7r + 7r _ mode, fo selected 



6.6 ±0.5 
2.45 ±0.31 
0.76 ±0.09 

6.0 ±0.6 
0.23 ±0.07 
0.94 ±0.09 
0.56 ±0.16 
1.83 ±0.24 a 
0.18 ±0.04 6 
0.17 ±0.07 c 
0.72 ±0.13 d 

0.75 ±0.09 
no entry 
0.060±0.014 
0.117±0.029 
0.068±0.024 e 



B 



V>(2S)^/ 



fa 



togram in Fig.[407a) is from the <p sidebands of Fig.[l"27c). 
and represents the dominant background in the (pir + n~ 
mode. The hatched histogram in Fig. l4"0T b') is from 
the X2K2tt° control region, and represents the dominant 
background in the </>7r°7r° mode. Subtracting these back- 
grounds, and subtracting a small remaining background 
using J/tp or tp(2S) sideband events, we find 181±19 
J/tp — > cf>ir + ir~ events, 45±9 J/tp — > <pir°ir° events, and 
19±6 tp{2S) — > 4>it + it~ events. We convert these to 
branching fractions and, after correcting for the modes 
other than <p -> K+K~, list them in Table [XlTTl All are 
consistent with current PDG values, of which the first 
two are dominated by our previous measurement. 

We do not observe any evidence for F(4260) decays to 
these modes, nor do we see a Y(4260) signal in any other 
mode studied here. 

Figures HIT a) (b) show the corresponding mass dis- 
tributions for 0/o (980) events, i.e., the subsets of the 
events in Figs. HDTa) and HDTb) with a di-pion mass in 
the range 0.85-1.10 GeV/c 2 . Signals at the J/tp mass 
are visible in both cases. From Fig. l4"TT b) we estimate 
16 ± 4 (f>f events in the 7r°7r° mode. However, 0/ o (98O) 
is not the dominant mode contributing to J/tp — >4>n + Tr~ 
decay. The open histogram of Fig. shows the 7r + 7r~ 



invariant mass distribution for events in the J/ip peak of 
Fig.HOta) (\m(K+K-ir+TT-)-m(J/iP)\ < 0.05 GeV/c 2 ); 
events in the J/tp sidebands (0.05 < \m(K + K ~7r + 7r~) — 
m(J/tp)\ < 0.1 GeV/c 2 ) are shown by the hatched his- 
togram. A two-peak structure is visible that is very sim- 
ilar to that studied by the BES Collaboration [31| and 
observed in Df — » 7r + 7r~7r + decay [32]. In both cases 
the 7r + 7r~ system is believed to couple to an ss system; 
both 7r + 7r~ distributions exhibit a clear /o(980) peak and 
a broad bump in the 1.3-1.5 GeV/c 2 region. The analy- 
sis of Refs. [U [32| shows that this bump is made up of 
/2(1270) and /o(1370) contributions; we denote this re- 
gion by f x . By selecting / (980) in the 0.85-1.10 GeV/c 2 
range and f x in the 1.1-1.5 GeV/c 2 range, shown by ver- 
tical lines in Fig. 221 and subtracting J/tp sideband back- 
ground we find 57 ±9 J /tp^-cpf {980) events and 61 ±10 
J/tp — > 4>f x events. 

Using Eq. (|10[) and dividing by the appropriate branch- 
ing fractions, we obtain the J/tp branching fractions 
listed in Table IXIIIl The measurements of Bji^^^fo in 
the 7r + 7r~ and ir°ir° decay modes of the fo are consistent 
with each other and with the PDG value, and combined 
they have roughly the same precision as given in Ref . [H . 

Note that, in contrast to </>(1680) — > (pirn decay, there 
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is no indication of a J ftp — > 0/o(6OO) decay mode. Only 
J /if) — > 0/ o (98O) is observed, as is true for the F(2175) 
state. 

We also observe 12 ±4 tjj(2S)—>^)fo, / ^> 7r+7r ~ events, 
which we convert to the branching fraction listed in Ta- 
bic IXIII) it is consistent with the value in Ref . Q , assum- 
ing B fo ^ v +„- = 2/3. 

The hatched histogram in Fig. [301a 1 ) shows the K + K~ 
invariant mass distribution, when the other kaon pair 
is in the (f> region, for the K + K~ K + K~ events in the 
J/t/j peak, selected by requiring \m(K + K~K + K~) — 
m{J/ip)\ < 0.05 GeV/c 2 . Subtracting sideband events we 
find 163 ± 19 events corresponding to J/ip — > <f>K + K~ 
decay. Using our normalization we obtain the branch- 
ing fraction listed in Table IXIII1 which agrees with that 
in Ref. Q but has better precision. In obtaining these 
values, we have used B(<p K + K~) = 0.489 [H], and 
assume equal rates for J/i/j — > (f>K + K~ and J/ip — > 
4>K°K°. 



X. SUMMARY 

We use the excellent charged-particle tracking, track 
identification, and photon detection of the EABAR 
detector to fully reconstruct events of the type 
e + e~ — > r ye + e~ — > ~fK + K~7r + ir~ , jK + K~ir°Tr° , and 
7_ftT + K~ K + K~ , where the 7 is radiated from the ini- 
tial state e + or e~ . Such events are equivalent to direct 
e + e~ annihilation at a cm. energy corresponding to the 
mass of the hadronic system. Consequently, we are able 
to use the full BABAR dataset to study annihilation into 
these three final states from their respective production 
thresholds up to 5 GeV cm. energy. The K + K~iy + tt~ , 
K + K~ir°ir° and K + K~ K + K~ measurements are con- 
sistent with, and supersede, our previous results 0- 

The systematic uncertainties on the e + e~ — > 
K+K-TT+7T-, K+K-TT 7r° and K+K'K+K- cross sec- 
tion values are 4%, 7% and 9%, respectively, for E c . m . < 
3 GeV, and increase, respectively, to 11%, 16% and 13% 
in the 3-5 GeV range. The values obtained are consider- 
ably more precise than previous measurements, and cover 
this low-energy range completely. As such they provide 
useful input to calculations of the hadronic corrections 
to the anomalous magnetic moment of the muon, and of 
the fine structure constant at the Z° mass. 

These final states exhibit complex resonant substruc- 
tures. For the K + K~iy + tt~ final state we measure 
the cross sections for the specific channels e + e~ — > 
K*(892) K-7T+, 07T+7T-, and cf>f , and, for the first 
time, for the e+er ->■ A 2 *(1430)°A-7r+ and e+e~ -> 
p(770)°K + K~ reactions. We also observe signals for the 
J\i(1270), Ai(1400), and / 2 (1270) resonances. It is diffi- 
cult to disentangle these contributions to the final state, 
and we make no attempt to do so in this paper. We note 
that the p° signal is consistent with being due entirely to 
Ki decays, and that while the total cross section is domi- 
nated by the A*(892)°A-7r + channels, only about 1% of 



the events correspond to the e+e~^ A*(892)°A*(892)° 
reaction. 

For the K + K~ -k°it final state we measure the cross 
section for e + e~ — > <fifa, and observe signals for the 
A*(892) ± and X^(1430) ± resonances. Again, the total 
cross section is dominated by the K *(892) + A~7r° chan- 
nel, but about 30% of events are produced in the e + e~— > 
A*(892)+A*(892)~ reaction. For the K+R-^tt final 
state we note that the cross section is roughly a factor 
of four smaller than that for K + K~ tt + h~ over most of 
the E c m . range, consistent with a factor of two isospin 
suppression of the 7r°7r° final state and another factor of 
two for the relative branching fractions of the neutral and 
charged K* to charged kaons. 

With the larger data sample of the present analysis, we 
perform a more detailed study of the e + e~ — > 0(lO2O)"7T7r 
reaction. The 7r + 7r~ and 7r°7r° invariant mass distribu- 
tions both show a clear /o (980) signal, and a broad struc- 
ture at lower mass interpreted as the /o(600). We ob- 
tain parameter values for these resonances. The 07r + 7r~ 
cross section measured in the K + K~it + ii~ final state 
shows a structure around 1.7 GeV and some structures 
above 2.0 GeV. The corresponding 07r°7r° cross section 
in the A + A _ 7r°7r° final state shows similar behavior. If 
the /o(980) is excluded from the di-pion mass distribu- 
tion, no structures above 2.0 GeV are seen. We fit the 
observed cross section with the VMD model assuming 
0(1680) -> </>/o(600) and 0(1680) ->■ 0/ o (98O) decay; the 
latter appears to be responsible for the threshold increase 
of the cross section at 2.0 GeV. Confirming our previ- 
ous study 0|, our data require an additional resonance 
at 2.175 GeV, which we call the V(2175), with decay to 
0/o (980), but not to </>/o(600). Further investigation re- 
veals consistent results for the K + K~ K + K~ final state, 
and clear V(2175) signals in the A+A~/ (980) channels, 
with /o(980)— ^ 7r + 7r~ and 7r°7r°. This structure can be 
interpreted as a strange partner (with c-quarks replaced 
by s-quarks) of the V(4260) [35| . which has the analo- 
gous decay mode J/i/jtt + tt~ , or perhaps as an ssss state 
that decays predominantly to </>/o- 

In the K + K~K + K~ mode we find e+e~^ 4>K + K~ 
to be the dominant channel. With the current data sam- 
ple we can say little about other K + K~ combinations. 
We observe an enhancement near threshold, consistent 
with the cj>fo channel and if these events are selected we 
have an indication of a V(2175) signal. Two other struc- 
tures in the K + K~ invariant mass spectrum are seen: 
the smaller could be an indication of the 0/o(137O) final 
state, and the larger of the (1525) mode. If events cor- 
responding to the 0/2(1525) final state are selected, the 
K + K~ K + K~ cross section shows a resonance- like struc- 
ture around 2.7 GeV, and a strong J/ip signal, which has 
been studied in detail by the BES Collaboration [3l|. In 
the K + K~ K + K~ cross section we observe a sharp peak 
at 2.3 GeV, which corresponds to the cj>K + K~ channel 
with the K + K~ invariant mass in the 1.06-1.2 GeV/c 2 
region. 

We also investigate charmonium decays into the stud- 
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icd final states and through corresponding intermediate 
channels, and measure the product of the electron width 
and the corresponding branching fraction. Some of the 
obtained J/ip branching fractions listed in Table IXIIII 
are as precise as, or more precise than, the current world 
averages, many of which were obtained in our previous 
study 0; the latter are superseded by our new results. 
We do not observe the K(4260) in any of the final states 
examined. 



ACKNOWLEGMENTS 

We are grateful for the extraordinary contributions of 
our PEP-II colleagues in achieving the excellent lumi- 
nosity and machine conditions that have made this work 
possible. The success of this project also relies criti- 
cally on the expertise and dedication of the computing 
organizations that support BABAR. The collaborating in- 
stitutions wish to thank SLAC for its support and the 
kind hospitality extended to them. This work is sup- 



ported by the US Department of Energy and National 
Science Foundation, the Natural Sciences and Engineer- 
ing Research Council (Canada), Institute of High Energy 
Physics (China), the Commissariat a l'Energie Atom- 
ique and Institut National de Physique Nucleaire et de 
Physique des Particules (France), the Bundesministerium 
fur Bildung und Forschung and Deutsche Forschungsge- 
meinschaft (Germany), the Istituto Nazionale di Fisica 
Nucleare (Italy), the Foundation for Fundamental Re- 
search on Matter (The Netherlands), the Research Coun- 
cil of Norway, the Ministry of Science and Technology 
of the Russian Federation, and the Particle Physics and 
Astronomy Research Council (United Kingdom). Indi- 
viduals have received support from jjjjjjj acknowledge- 
ments. tex CONACyT (Mexico), the Marie-Curie Intra 
European Fellowship program (European Union), the A. 
P. Sloan Foundation, the Research Corporation, and the 
Alexander von Humboldt Foundation. ======= the 

Marie-Curie IEF program (European Union) and the A. 
P. Sloan Foundation (USA). 

ilLLLLL 1-24 



V. N. Baier and V. S. Fadin, Phys. Lett. B27, 223 (1968). 

A. B. Arbuzov et al, J. High Energy Phys. 9812, 009 [19] 
(1998). 

S. Binner, J.H. Kiihn and K. Melnikov, Phys. Lett. 
B459, 279 (1999). [20 
M. Benayoun et al, Mod. Phys. Lett. A14, 2605 (1999). [21 
C. K. Nakamura et al. (Particle Data Group), Journal of 
Physics G37, 075021 (2010). [22 
See for example: S. Gomez-Avila, M. Napsuciale, E. 
Oset, Phys. Rev. D79, 034018,2009; M. Napsuciale, E. [23 
Oset, K. Sasaki, C.A. Vaquera-Araujo, Phys. Rev. D76, 
074012 (2007). [24 

B. Aubert et al. (BABAR Collaboration), Phys. Rev. D76, 
012008 (2007). [25 
M. Ablikim et al. (BES Collaboration), Phys. Rev. Lett. [26 
100, 102003 (2008). 

CP. Shen et al. (Belle Collaboration), Phys. Rev. D80, [27 
031101(R) (2009). 
M. Davier, S. Eidelman, A. Hocker, Z. Zhang, Eur. Phys. [28 
J. C31, 503 (2003). 

B. Aubert et al. (BABAR Collaboration), Phys. Rev. D69, [29 
011103 (2004). [30 
B. Aubert et al. (BABAR Collaboration), Phys. Rev. D70, 
072004 (2004). [31 
B. Aubert et al. (BABAR Collaboration), Phys. Rev. D71, 
052001 (2005). [32 
B. Aubert et al. (BABAR Collaboration), Phys. Rev. D73, 
052003 (2006). [33 
B. Aubert et al. (BABAR Collaboration), Phys. Rev. D76, 
092005 (2007). [34 
B. Aubert et al. (BABAR Collaboration), Phys. Rev. D77, 
092002 (2008). [35 
B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 
103, 231801 (2009). 

A. Cordier et al. (DM1 Collaboration), Phys. Lett. B110, 



335 (1982). 

B. Aubert et al. (BABAR Collaboration), Nucl. Instrum. 
and Meth. A479, 1 (2002); W. Menges, IEEE Nucl. Sci. 
Simp. Conf. Rec. 5, 1470 (2006). 

H. Czyz and J. H. Kiihn, Eur. Phys. J. C18, 497 (2001). 
A. B. Arbuzov et al, J. High Energy Phys. 9710, 001 
(1997). 

M. Caffo, H. Czyz, E. Remiddi, Nuovo Cim. A110, 515 
(1997); Phys. Lett. B327, 369 (1994). 
E. Barberio, B. van Eijk and Z. Was, Comput. Phys. 
Commun. 66, 115 (1991). 

S. Agostinelli et al. (GEANT4 Collaboration), Nucl. In- 

str. Methods Phys. Res., Sect. A 506, 250 (2003). 

T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 

S. Jadach and Z. Was, Comput. Phys. Commun. 85, 453 

(1995). 

N.N. Achasov and A. A. Kozhevnikov, Phys. Rev. D55, 
2663 (1997). 

A. Aloisio et al. (KLOE Collaboration), Phys. Lett. 
B537, 21 (2002). 

S. Pacetti, Eur. Phys. J. A31, 665 (2007). 

T. A. Armstrong et al. (WA76 Collaboration), Z. Phys. 

C52, 389 (1991). 

M. Ablikim et al. (BES Collaboration), Phys. Lett. 
B607, 243 (2005). 

B. Aubert et al (BABAR Collaboration), Phys. Rev. D79, 
032003 (2009). 

M. Ablikim et al. (BES Collaboration), Phys. Lett. 
B598, 149 (2004). 
M. Ablikim et al 
B698, 183 (2011). 

B. Aubert et al. (BABAR Collaboration 
95, 142001 (2005). 



(BES Collaboration), Phys. Lett. 



Phys. Rev. Lett. 



